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SECTION I 


INTRODUCTION 

The current state of the design of physical systems is 
in some respects as much an art as a science. Of the many 
considerations ' involved in systems design, perhaps those 
related to an evaluation of the physical performance best 
lend themselves to formal automation. This report describes 
the results of an effort to place this aspect of real-world 
systems design on an analytical basis. In particular, a 
computational capability, PERFORM, has been developed for the 
evaluation of the limiting physical performance of systems 
subject to transient disturbances. 

PERFORM provides the characteristics of the theoreti- 
cally best, that is, the limiting , design concept according 
to response criteria. As a consequence, for certain mechani- 
cal, structural, and control systems PERFORM makes it pos- 
sible to approach a design directly from the design criteria 
with no a priori commitment to a particular design concept. 

With this limiting performance capability, the designer can 
determine the feasibility of his proposed design on the basis 
of the specifications alone; moreover, he can monitor and 
measure his success during the design process itself. Here- 
tofore, without the characteristics of the limiting design, 
the evaluation of proposed designs could be made only by per- 
forming a multitude of analyses for each candidate design. 

This capability applies to systems with transient load- 
ing in which response criteria are of concern. These criteria 
are usually expressed in terms of constraints formed from 
displacements, stresses, forces, velocities, or accelerations 
that cannot exceed certain prescribed values. In addition, 
objective functions of peak values of other responses are to 
be minimized. PERFORM provides the time-optimal character- 
istics of portions of a system such that these constraints 
are satisfied and objective functions are minimized or maximized . 

An outline of the capabilities of PERFORM is given in 
Section II. 

Section III of this report describes the components and 
design of PERFORM. Also, the procedure for using PERFORM is 
outlined. 

Section IV contains applications of PERFORM, including 
some numerical results. Among the problems discussed are rail 
vehicle suspension systems, train impact, aircraft ride con- 
trol, launch vehicle control, missile/silo isolation, and a 



nuclear reactor control system. 

The Users Guide for PERFORM is given in Appendix I. The 
details of the technical formulation underlying PERFORM are 
in Appendix II. Appendix III contains the programming aspects 
of the system, including programming documentation; The final 
Appendix contains listings of programs. 



SECTION II 


CAPABILITIES OF PERFORM 

PERFORM is a computer system that can be used to deter- 
mine the limiting performance characteristics of a dynamic 
system subject to transient loading. ,Ref. 1 considers in 
detail the concept of limiting performance and its application 
to shock isolation systems. The example problems of Section 
IV provide an indication of the range of applications of PER- 
FORM. 

The dynamic system can be described and input to PERFORM 
using the first or second order equations 

s = As + Bu + Df ^ 

•• • 

Mq + Cq + Kq + Uu = Ff , 

in which u is a vector of time varying functions, called con- 
trol or isolator forces, that have replaced portions of the 
dynamic system. A, B, D, M, C, K, U, F are coefficient mat- 
tices. s and q are vectors of response varables, e.g., dis- 
placements, stresses, accelerations. f^ is a forcing function 
vector where the subscript k designates K the k t set of forcing 
or loading functions. This allows the system to encounter 
alternative sets of disturbances which might occur with equal 
probability. 

The acceptable equations of motion appear to be linear. 

In fact, however, they are_"quasilinear " since those portions 
of the system replaced by u can be linear, nonlinear, active, 
or passive. The remainder of the system must be linear as 
must the overall kinematics. 

The user must place his equations in one of the forms 
of the above equations. The non-zero elements of the mat- 
rices A, B, D or M, C, K, U, F are then entered as inputs. 

This is accomplished by identifying the matrix, e.g., M MAT- 
RIX, and then specifying an element and its value, e.g., i, 
j, and M^ j . Elements not entered are assumed to be zero. 

PERFORM finds the characteristics, including u and trade- 
offs between optimal response variables, of the dynamic system 
such that bounds on some of the response variables s or q or 
control forces u are not violated while the maximum (or mini- 
mum) in time of other elements of s or q are minimized (or 
maximized) . 
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Regardless of the form (first or second order) used to 
describe the equations of motion, the formats for the objec- 
tive function and constraints are the same. In the case of 
the second order equations, a state variable vector § is esta- 
blished as 



Any linear combination of state, variables, derivatives of 
state variables, or control forces can be used as an objective 
function. In the case of the system described by second, order 
equations, these become linear combinations of accelerations, 
velocities, displacements, and control forces. The objective 
function is input to PERFORM in the form of 

PXl s + PX2u +PX3? k 


where PXl , PX2 and PX3 are coefficient matrices. If more; than 
one row of the matrices of this equation contains non-zero . 
elements, then the peak values in time of the vectors result- 
ing from the meaningful rows are to be compared. PERFORM min- 
imizes (maximizes) the maximum (minimum) of the peak values. 

Constraints may be placed on state variables, derivatives 
of state variables, and control forces. The general form, 
which is again linear, is - . 

YL £ Yl s + Y2u + Y3f k £ YU 

where Yl, Y2, Y3_ are coefficient matrices and YL, YU are lower 
and upper bound vectors. Constraints can be imposed at every 
time of the response or at specific times. 

In summary, PERFORM accepts system equations of the form 
given above. _For prescribed initial conditions,. PERFORM then 
-computes the u vector such that the max | PXl s + PX2u + PX3 f ^ [ 

is minimized ( or min | | is maximized) while the above con- 

straints are satisfied. Any linear combination of s, u, and 
f k can be tabulated or plotted versus time. A tradeoff curve 
between the maximum objective function and any particular con- 
straint can be generated by varying the bounds on that con- 
straint. 
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1 ' -- ' SECTION III 

'Zj . j t 

' PERFORM SYSTEM; DESCRIPTION ■ ■ • 

To do the limiting performance problem calculations a 
computer system named PERFORM has been developed. Because 
several separate programs are used, it is referred to as a 
system rather than a program... 

In this section the system is described with increasing 
detail. First, the system is presented from the "systems 
analyst'!' point of view; the programs, data sets, and inter- 
relations, thereof are described and graphically presented with 
a system.', 1 flowchart'.' Next the functions, subroutine calls, 
logie7,~ O.tc'. of each of the system programs are described. Pro- 
gram' listings and discussions are in Appendix III. 


System Design 


PERFORM , is a system of several programs, the functions, 
inputs, ;and outputs: of which are outlined here. Influencing 
considerations, of., the development and operation environment 
are-- r a r lso discussed. 

— J'.O MtfCV.T; ■; 


cf>^l .V. 


astern Flow Description 


Figure 1 is the PERFORM system flowchart. The three 
separate system programs are represented by rectangular boxes; 
these programs are: 


. PREPROC (Pre-processor) 

. LP Solver (linear programming computer 
program) 

. PSTPROC (Post-processor) 
a. PREPROC ■ 


7 ’ Using the mathematical model of the dyna- 

mic .-^system, the pre-processor program computes and punches 
the ''linear programming problem from the PERFORM problem sped-*- 
fications described in Appendix I. There are three outputs of 
this’' program: 

' . PSTPROC Report Specification. This card 

file contains information necessary for 
■ , PSTPROC (to which this is an input file) 

•• to determine and print the problem solu- 

tion. 


Linear Programming Problem Input. 
This card output is the linear pro- 
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SPECIFICATIONS 

DATA 



Fig* I PERFORM System Flowchart 
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gramming problem input in MPS/360 format, 
which is accepted by most major LP solvers. 

PREPROC Report. This report is a printed 
record of the above files. 

User options can control the contents of 
this report. The curved-bottom box is the 
standard flowcharting symbol for printed 
reports 

The data transmission links are restricted somewhat in 
both the University of Virginia development environment and the 
NASA/Langley operating environment. These constraints are dis- 
cussed in the section entitled "Environment Considerations." 

PREPROC is written entirely in FORTRAN and has been tested 
on the University of Virginia Burroughs B5500 ,CDC6400 , and the 
NASA/Langley CDC 6600. 

b. LP Solver 

This refers to any software package that 
can solve large linear programming problems. Examples of such 
capabilities are IBM MPS/360, CDC OPTIMA, and CDC OPHELIE . 

All of these programs accept MPS/360 input format. Since MPS/ 
360 was used in the developmental stages of PERFORM, the LP 
Solver in the subsequent discussions of this report will be 
referred to as MPS/360. Anywhere MPS/360 is mentioned, other 
software packages, such as OPTIMA, could also be used. 

MPS/360, Mathematical Programming System/360, is an open 
ended IBM software package capable of efficiently solving 
extremely large linear programming problems. The latest ver- 
sions can handle problems of up to 16,000 rows and essentially 
unlimited columns. Several workfiles used by MPS/360 are omit- 
ted for clarity on the system flowchart of Fig. 1. 

The input to MPS/360 is the output of "linear programming 
input" of PREPROC. 

C. PSTPROC 


The post-processor program produces the 
final report of the problem solution. Two inputs are; neces- 
sary for PSTPROC: the linear programming problem solution 

from LP Solver and the PSTPROC report specif ications from 
PREPROC. This program is written entirely in FORTRAN. 

2. Environment Considerations 

In Fig. 1 the programs are labeled "local" and "external" 
computer and a communications link is shown to and from the LP 
Solver step. This clearly awkward arrangement is used because 


7 



of the limitations of available computing facilities. The 
Langley Research Center CDC 6600 computing facility has no 
appropriate linear programming software. At Langley the PER- 
FORM system cannot be run entirely, locally, thus the need for 
an external computer. The system. developed is a synthesis of 
these above constraints and the desire for maximum utiliza- 
tion of local facilities. - ; - 

Punched cards are indicated '.as the storage medium for 
program data interchange. Cards, were convenient for trans- 
mission to the external computer during development. Replace- 
ment with tape would be desirable and easily accomplished for 
large problems with unmanagable card files. ]. 

t 

The decision to use MPS/360 was based on economy ahd 
availability. For development, the Virginia- Polytechnic 
Institute and State University IBM/360 50-65 multiprocessor 
system, with an IBM 2780 telecommunications terminal was used. 
Recently, OPTIMA on the University of Virginia CD C 64001 has 
been used. i 


Use of the system is much simpler if a linear program- 
ming package is available locally. The data can be passed 
between programs on a disk, and the user does not have t‘o be con 
cerned with this data. One can simply put his problem speci- 
fication data with some prepared control cards and receive 
the problem solution with one run. * 

B. Program Descriptions 

Each of the three programs of the system is further 
described in this section. The purpose of the description 
is to relate the program coding to the limiting performance 
problem model but without involving computational details. 

A detailed listing with more discussion is given in Appendix 
III. 

1. PREPROC 

Figure 2 is a flowchart of the main section of 
PREPROC. This section contains all input, calls to computa- 
tional subroutines, and varying bound calculations. 

Some notes about the flowchart are in order. Predefined 
blocks, the blocks with outward pointed sides and parallel top 
and bottom, represent subroutine calls; the name in the block 
is that of the subroutine. The trapezoidal-shaped box enclos- 
ing most of the third page and beneath the label "For each 
varying bound" is unconventional. It signifies that the se- 
quence of instructions in the block is to be executed for the 

* This is currently possible on an automated version of PERFORM at the 
University of Virginia Computer Science Center, both for on-campus users 
and for users accessing the computer from a remote, compatible computer 
terminal. 8 
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> Fig. 2. PREPROC Main Section 





Fig. 2 PREPROC Main Section (continued) 
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pig. 2 PREPROC Main Section (concluded) 
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value of each varying bound, i.e., varying row of the YL and 
YU vectors. This flowchart is not an exact representation of 
the actual program; from it one can tell what the program does 

but not necessarily how it is done. 

>. 

Because all files except those for printed output in the 
PERFORM system are either cards or in card format, the word 
"record" used to describe file elements can be interpreted as 
"card" in this discussion. The records of the problem speci- 
fications file are of- two types: data and indicator. Data 

records contain values of program variables. By preceding 
data records, indicator records associate values with varia- 
bles; indicator cards also contain problem processing speci- 
fications. 

Before each section of data records of values for pro- 
gram variables there is an indicator card specifying the pro- 
gram variable array. The last card of a section contains a 
1' in the key field; previous cards of the section have a zero 
in this field. Following the last card of a section the pro- 
gram presumes an indicator card . The program contains no 
error detection facilities so it behooves the user to double 
check his input. Details of these records are found in Appen- 
dix I. The first thing done by the program is to read three 
problem specifications records; in Appendix I these three 
records are numbered cards 1, 2, and 3. Next, most program 
variables are initialized to zero. 

The fourth block reads an indicator card. If the card 
is a STOP indicator card the program stops. If a RESTART in- 
dicator, control returns to the program beginning and a new 
problem may be processed. A FINISH INPUT DATA record starts 
the problem processing. Appropriate processing for other in- 
dicators is performed - this includes such things as array 
entry, array printing, problem modification, etc. 

If the dynamic system equations are second order, sub- 
routine COVERT is called to reduce the second order equations 
to first order. The order reduction method is described in 
Appendix II. 

Subroutine MEXAP computes exponential functions of mat- 
rices. These are array P and array AEI , both defined in 
Appendix II. 

The next routine , BRIDGE , computes arrays for matrices R, 

T, and vectors Cl, C2, and C as described in Appendix II. 

Routine ELEMTP is called to compute those portions of 
the H and G arrays that are determined by the objective func- 
tion. These portions are shown in Table 1 of Appendix II. 
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The other portions of these arrays, determined by the prob- 
lem constraints, as shown in Tables 2A and 2B of Appendix 
II, are calculated in subroutine ELEMTQ. Subroutines ELEMTR, 
ELEMTS, and ELEMTT are called by ELEMTQ for calculations but are 
not shown in Fig. 2. With H and G arrays computed, the next 
step is to output the results. 

Subroutine LPDATA is called to produce the MPS/360 
linear programming problem. Other subroutines, not shown, 
are called to punch different sections of the MPS input:. 

CODATA for columns, RODATA for rows, RHDATA for the right hand 
sides, and BODATA on the bounds. 

The POST routine punches the PSTPROC report specifica- 
tions. Although this file is shown separately from the MPS/360 
linear , programming file the user must separate the two from 
the single card punch output. 

To provide a printed problem record subroutine PRNT prints 
the H and G arrays. If a printout of the calculated arrays is 
desired, as indicated by a nonzero value of variable INT . in / 
the VERIFY section, it is produced by subroutine VERFY. The 
problem specification is produced when the VERFY section is . 
processed and variable INP is nonzero. 

Any number of the lower or upper bounds of the rows of 

the matrix expression YU<_ lYL may be varied. 

Appendix I contains details of how this is specified. 

For each value of each varying bound the linear program- 
ming problem is computed, punched, and printed by the sub- 
routine sequence ELEMTP , ELEMTQ, LPDATA, and PRNT. If a print- 
out of the problem arrays is desired as indicated by a nonzero 
value of variable INT, it is produced by subroutine VERFY. 

When the loop has been performed for each value of each vary- 
ing bound, control returns to the indicator point for modifica- 
tions or other instructions. 

There are several additional array manipulation subrou- 
tines in PREPROC not mentioned here. These are described in 
Appendix III A. 

2. MPS/360 


Mathematical Programming System/360 is an open ended 
IBM software package for solving optimization problems. MPS 
is not a program, but a set of procedures with linking facili- 
ties. The procedures can do such things as read data, set up 
files, write solutions, and, of course, solve linear program- 
ming problems. The procedure linking facility is known as 
the control program. It consists of procedure call statements 
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in addition to ordinary programming facilities such as GOTO, 

IF, MOVE statements, etc. This control program is compiled 
and executed in separate OS job steps. 

There are two outputs from MPS: the printed linear pro- 

gramming problem report and the punched problem solution. 

If the punched linear programming problem output of PRE- 
PROC is interpreted, a similarity to the PREPROC input can be 
observed. Nonzero program variables are entered on cards with 
array row and column numbers and follow, section identifying 
cards. Unlike PREPROC the alpha characters ROW, COL precede 
the left justified numbers. The characters are necessary be- 
cause MPS identifies rows and columns of arrays with FORTRAN- 
like variables which must begin with an alphabetical character. 
ROW and COL serve to identify columns of the linear program- 
ming tableau; RHS identifies the right hand side vector' and 
BDS identifies the variable bounds specifications. • 

.... o 

The control program of MPS/360 for PERFORM is rather sim- 
ple; the data is read in, the problem solved, and the solution 
written and punched. A listing and discussion of the control 
program is in Appendix IIIB. 

Each. linear programming problem produced for each value 
of each varying bound requires a separate MPS/360 run. 

3. PSTPROC 


Program PSTPROC writes the solution to the limiting 
performance problem using the MPS/360 linear programming pro- 
blem solution and the PSTPROC report data from PREPROC . 

Figure 3 is a functional flowchart of PSTPROC. First the 
PSTPROC report data, punched by subroutine. POST in program 
PREPROC, is read. Next the first linear programming problem 
solution is read. Again it should be pointed out that no error 
detection facilities are in the program. If a card is drop- 
ped from the PSTPROC report specifications a linear program- 
ming solution card will take its place. 

For each row of the tabulation expression, i.e., the Q1 , 
Q2, and expression, the block "inside the box is executed. 

For all time intervals the value of the expression is calcu- 
lated and printed. If a graphical presentation of these cal- 
culations is desired, as specified by a nonzero value ITR(I) 
where I is the row number in Q^, Q2_, Q3; -matrix expression, 
subroutine TRAJ is called. 

After all the tabulation expressions have been procured 
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Fig. 3 PSTPROC Main Section, Functional, Fow Chart 


15 



the next linear programming solution' is read and the calcula- 
tion loop is repeated. If all solutions have been processed 
variable ITD is checked. . ’ 

If ITD is nonzero, a trade-off diagram is printed by 
subroutine TRADE. 

After the trade-off diagram is produced processing stops. 
Note that PSTPROC handles only one problem -a ; run; different 
problems, with different PSTPROG report decks , require sepa- 
rate runs. . .. ' , ;!• ' . 
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• • . SECTION IV . . i< ’ 

•« . "< t t. • I iJ i x . i '• 

APPLICATIONS 

A. Single-Degree-of-Freedom Shock Isolation ■ System 

We begin the applications of PERFORM by considering a 
single-degree-of -rfreedom (SDF) dynamic system subject to . tran- 
sient^ loading. Since this system has been thoroughly explored 
elsewhere (e.g. , Ref. 1) , it is included here only to demon- 
strate the use of PERFORM in the study of the limiting perfor- 
mance of a familiar system. 

Consider the SDF system together with the acceleration 
shock shown in Fig. 4. The equation of motion is given by 

mz + u = 0 (1) 

or, for a unit mass, 

z + u = 0 (2) 

with initial conditions z(0) = z(0) = 0 
The kinematic relations are 

z = x + y (3a) 

z = x + y (3b) 

z=x+y = x + f (3c) 

Substituting Eq(3c) into Eq(2) gives 
x + u = -y = -f 
x(0) = y (0) =0 
x (0 ) = y (0) =0 

Suppose we wish to find the lowest possible peak acceler- 
ation of the mass (z) if the relative displacement between the 
mass and the base is constrained. The function u(t) is the 
system controller or the isolator function. Thus, we wish to 
find u(t) such that the maximum acceleration of the mass 

<J> = max | z | = max | u | 

is minimized while the relative displacement satisfies the con- 
straint 

. |x| < A 
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We seek the trade-off between minimum <p and a range of values 
of A. The design value of this information is discussed in 
Ref . 1 . . 

Specifically, we choose to compute min <{> for A = 0.01, 
0.02, and 0.03. 

The formulation of the problem in the PERFORM format 
follows. The notation is defined in Appendix I. 

1. Equations of Motion (x + u = -f) 

• • • 

Mq + C q + Kq + Uu = Ff ^ 

where 

M = [1] 

C = 

K = t0] 

U = [1] 

F = t"l] 

q = M 

‘ q = [x] 

5 = [x] 

u = [u] 

? k = [y] = [f] 

2. Forcing Function (Fig. 4b) 

f (l,j,l) = 10.0 for j = 1 to 5 

f(l,j,l) = 0 for j = 6 to 9 

where the first index 1 designates that this is the first 
(and the only) set of forcing functions, j = discrete time, the 
last index 1 indicates the first component in forcing function 
set No. 1. 

3. Objective Function (4)) 

PX1 s + PX2 u + PX3 f 
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where 



PX1 = 

to 

s = 

X 


X 

PX2 = 

tl] 

PX3 = 

[0] 


4. Constraints (|x|£A) 

YL < Yls' +' Y2u + Y3f k <_ YU 


% 


where 

Y1 = [0 1] 

Y2 = [0] 

xi = to] 

Since the constraint is for all time and of (<_ <) type we 

set 


MSP (1,1) = 0 
MSP (1,2) = 0 

To obtain a trade-off diagram the value of A is varied. 

Thus 

YU = A q1 + nAA x ' 

YL = A 02 + nAA 2 

where 


ii 

i — i 

o 

< 

SVA1 =0.1 

li 

CM 

O 

C 

SVA2 = -0.1 

AA^ = 

0.1 

t> 

> 

to 

II 

-0.1 

n = 

1 to NIC 
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5. Output 


The time trajectory of u is obtained .using the out- 
put matrix 


Qls + Q2u + Q3f k 

Q1 = [00] ■, j 

02 = [ 1 ] 

03 = [03 : , . 

6. Other Data Pertinent to This Problem 


We set 

NOE = 
MODE , = 
NDF = 

NU =r 
NF = 
NSETS = 
II = 

■ « • » 

TM = 

NOB = 
NOC = 
ISP = 
NOT = 

and 

ISET(l) = 
The whole input 


2 • 

1 

1 , 

1 

1 

* 1 . i 

1 

10 

• 1 ' ■ 

0.02 

1 ; 

1 . 

0 

1 (here we assume only one output is 
required) ■ • . . 

1 

{ .* 

deck is shown in Fig. 5.; 
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7. Results 


i i i 

Figures 6a , through 6g show the output of PERFORM 
resulting from the input deck of Fig. 5. These are the time 
trajectories of the control forces ''(u) (Figs. 6a to 6f) and 
a diagram (Fig. 6g) indicating the performance trade-off 
between minimum <J> the objective function, and A, the bound 
of the constraint. 



_ — TBs. Q1*S*-Q2*U+_Q3*E- 
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. Qla 0. 00 0.00- 

02 = 1.00 

03s 0. 00 — 


TIHE INTERVAL 

TB 
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. 020 
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-. 02 0 _T 0_ 

. 040 

-ft. 3S3331 
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. 060 
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— ... . 080 - 
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.140 
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. . 140-TO. 

.160 
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.160 TO 

. 180 
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Fig. 6a. Trajectory of u(t) for |x|<0,01. 
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0.000 TO .020 

02 0— T 0 .040 

. 040 TO .060 

.060 TO... .080 

.080 TO .100 

.100 TO .120 
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Fig. 6b. Trajectory of u(t) for |x|<0.02. 
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1- TBs-01 *S* - Q2*U*-03*£. 

- WHERE 

01 = 0.00 0.00 

02 = 1.00 

— 03= 0.00 ...... — — 
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.000 TO 

.020 

-6.250C00 

_020 TO 

.040 

--6. 250000 

.040 TO 

.060 
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-6.20000 

.aJ.4Q TO 

.160 

... -6. '’0000 
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Fig. 6e. Trajectory of u(t) for | x | <0 .03 . 
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Fig. 6f. Plot of uCt) for |xj<0.03. 
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"'-OBJECTIVE FUNCTION ■ 8.333333 BOUND « .010100- 

OBJECTIVE FUNCTION * T.1NB05T BOUND ». .010000 

- OBJECTIVE FUNCTION ■ 6.Z50000 BCUNO » .030000 


Fig. 6g. Trade-off Diagram for SDF System 
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53 . Multi-Degree-of-Freedom Shock Isolation Systems 

To use PERFORM in evaluating the limiting performance 
characteristics of a shock isolation system the equations of 
motion of the system, after the control forces are introduced, 
must be linear. However, many dynamic systems are described 
by nonlinear equations and we are invariably faced with the 
problem of approximating the systems with various assumptions 
and yet not losing sight of the actual system. In this sec- 
tion, in addition to treating some multi-degree-of-freedom • 
shock isolation systems, we will consider some of the reason- 
able assumptions to be made for thes linearization of the equa- 
tions of motion for some dynamic systems . of nonlinear nature. 

1 . Klein’s Three-Degree-of-Freedom System 
a . Problem Description 

In Ref. 2 a three-degree-of-f reedom system as 
depicted in Fdg. 1 was considered. > The model to be analyzed 
has the z^ , z 2 plane as a plane of symmetry. It is assumed 

that the mass is a rigid body; that principal axes of inertia 
through the mass center parallel the edges of the body; and 
that the motion of the mass will hot affect the motion of the 
base. In addition, it is assumed that the points of appli- 
cation of the shock isolator forces u^, u 2 , u^ and u^ remain., 

fixed at the corners of the supported mass and retain their 
horizontal and vertical directions despite the motions of the 
base and mass. The kinematic conditions of the system are 
shown in Fig. 7, to which the following definitions apply: 

u. = the shock isolators or control forces which are 
to be optimized, 

y. = the two orthogonal input motions of the base referred 

to a fixed reference frame, 

z . = the two orthogonal motions of the center of mass 

relative to a fixed reference frame, 

0 = rigid body rotation of the isolated mass about the 

center of mass. 

The equations of motion for the system shown in Fig. 7 become: 

mZf - u^ + u^ = 0 

{4) 

mz' 2 - u 2 - u 3 = 0 
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Fig. 7 Active Isolation System Used In Ref. ! 
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10 - au 1 cos 0 + bu 2 cosG du 3 cos t 0, + au 4 cos 0 = 0 (4c) 

■ • ' . . . i. a:: ’ : . i "• . •: ■ 

where m is the isolated mass and I is the moment of" inertia' of 
the isolated mass about its center of mass. 

It is further assumed in . Ref. 2 that 0 is small , therefore 
allowing the usual small angle approximations to be made. Using 
the kinematic relations - . .. • . 


2 1 

= X 1 

+ 

y l' 

2 1 

= *1 

+ 

y l 

Z 2 

= 

+ 

y 2' 

5 2 

= *2 

+ 

y 2 


Eq. (4) becomes 

mXf - Uj. + u 4 = -my 1 

mx 2 - u 2 - u 3 = -my 2 ( 5 ) 

10 - au^ + bu 2 - du 3 + au^ =0 
2 

Letting I = mp , where p is the radius of gyration, we find 
for unit mass 

*i - u i + u 4 = ~h . 

x 2 ~ U 2 u 3 = -y 2 (6) 

0 - au. + bu 0 - du,. + au . = 0 ' . . 

p 2 r r r 

The above equations define the behavior of the isolated mass as 
a function of the applied accelerations y.^ and y 2 and the control 
forces ui, U 2 # U 3 , and U4. . , 

For this linearized system the optimization problem con- 
sidered in Ref. 2 is: Find the control functions u^t), u? (t) , 

Uq ( t ) and u^(t) so as to minimize some function of the maximum 
absolute values of the response acceleration, f(|z]J max 
I z 2 lmax, 1 9 Imax) when the base is subjected to input y. and/or 
y and the mass is constrained to move within the rattlespaces 
defined by and t 2 as shown in Fig. 7. 
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Since l : inear programming is selected as a solution tech- 
nique, the objective function and constraints must be expres- 
sible as, linear combinations of control forces. The objective 
function chosen was 

. . $ = max {max jw-^(t)|, max | W 2 ( t ) | , max jw^ (t) | } (7) 

-■ ■ 1 : i . . . • . 

* ; • I 

where w^ = (u^ - u^) = 

w 2 = “2 (u 2 + u 3> = a l S 2 ‘ (8) 



with a., a 2 as weighting constants. 

The constraints are 

l . » 

I X 1 + a0 ll e i 

|x^ - eGl^e^ 

l x 2 “ b0 ll e 2 O) 

|x 2 + d0|j<e 2 

The forcing function considered was 

y (t) = t 2 e _t 1 (10) 

This gives an acceleration 

, " ’ y (tj = e -t (t 2 - 4t + 2) (11) 

To compare the use of PERFORM with results available in 
Ref. 2, the forcing function is chosen to act in the y, direc- 
tion only. Other pertinent data include: 

No. of time intervals =40 

' • , . , #1 - 1.0 

a 2 = 1 '° 
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e 2 = 0.0 

a = 0.0 

b = 25.0 
d - 25.0 
e = 0.0 

b, Formulation in. PERFORM Format 


The Ref. 2 example outlined above will now be 
reformulated in accordance with the PERFORM format. Since the 
forcing function acts in the y, direction, u ' and u, will be 
inactive, as will w_. Hence only two functions will be includ- 
ed in the objective^ function. Due to e 2 =0, only the con- 
straints on will be taken into agcount. 

Equations of Motion 


Using the definitions for w,, w„ , and w~ as given in Eq. 
(8) and the values of the constraints as defined above, the 
equations of motion can be written as 

•• •• 

X 1 - w i = -Yj. 

x 2 - w 2 = 0 (12) 

0 - w^ = 0 

Eq. (12) may be rewritten in the PERFORM format as 

Mq + Cq + Kq + Uu = Ff R (13) 

where the vectors are defined as 





r 




"*r 


o 

X 1 

, 

X 1 

■ 


. 

• 



q = 

X 2 

q = 

X 2 

q = 

X 2 


0 


0 . 


0 


- 


- “ 


r. **■ 
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and the coefficient matrices as. 


10 0 
0 1 0 
0 0 1 

0 0 0 
C = | 0 0 0 

! 0 0 0 




0 0 0 
K = 0 0 0 

0 0 0 




• l 

The objective function given by Eq. (7) is to be minimiz- 
ed. In the PERFORM input format the objective function becomes 

PXls + PX2u + PX3f, (14) 

— — k 

where 
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( 15 ) 


Y1 = [0 0 0 1 0 0] 


Y2 = [0 0 0] . 

Y3 = [0] 

The input deck is shown in Fig. 8. 

c . Comparison of Results 

The results obtained by using a linear programming opti- 
mization procedure for this case given in Ref. 2 are dupli- 
cated in Table 1. The results of the same problem solved by 
PERFORM are given in Table 2a. The values of the objective 
function from these two solutions are very close to each 
other. The tabulated output from PERFORM is the wj in Ref. 2. 
Using the values of u^ and u 4 from Table 1, the w^ calculated 
from Eqs. (8) are about the same as those given by PERFORM for 
the first five time intervals. A comparison of the w^ values 
for the first 10 time intervals is given in Table 2b. Since 
the solution of the optimization problem is not unique after 
the maximum value of the objective function has been reached, 
it is not surprising that Table 1 and Table 2a do not agree 
consistently at later times. 


d. Modified System 

Based on certain assumptions, the above system can be 
modified to include large rotations. Additional isolator 
forces u and u R can be introduced for rotational motion of 
the isolated mals^ and rotational contributions of horizontal 

and vertical isolator forces can be neglected (i.e., omission 
of all but 0 in the third of Eqs. (6)). • The equations 
of motion for the system shown in Fig. 9 are written as 



x. 

- u + 

u , = 

-y 


1 

1 

4 

*1 


X„ 

- u - 

u = 

-y • (16) 


2 

2 

3 

2 2 


9 

+ u a + 

U 8 = 

0 . 

Eqs. (16) are 
sarily small. 

applicable 

to the 

case 

where 0 is not heces- 


2. Missile-Silo Isolation 


Next we consider another problem which is quite simi 
lar to the one discussed in Ref. 2 but has a different arrange 
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Fig. 
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. Input Deck for Multi-Degree-of-Freedom System 
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Fig. 8. Input Deck for Multi-degree-of-Freedom System (concluded) 
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Table 1. Results of Multi-Degree-of-Freedom System given in Ref. 2. 
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* 01*S« 

02 *U* 
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Table 2* a Results for Multi— Degree-of— Freedom System 
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i 

W 

u 4 (ti) 

wCt^ » u x -u 4 

u(t A ) 

1 

0.2948257 

-0.2948257 

0.5896514 

0.589653 

2. 

0.2848257 

-0.2948257 

0.5896514 

0.583650 

3. 

0.2948257 

-0.2948257 

0.5896514 

0.589650 

4. 

0.2948,257 

-0.2948257 

0.5896514 

0.589650 

5. 

0.2948257 

-0.2948257 

0.5896514 

0.589650 

6. 

0.0 

o' 

• 

o 

o 

• 

o 

0.117460 

7. 

-0.2948257 

0.2948257 

-0.5896514 

-0.589650 

8. 

-0.172954 

0.1727954 

-0.34559080 

-0.589650 

9. 

-0.2948257 

0.2948257 

-0.5896514 

-0.589650 

0 . 

0.0 

0.0 

0.0 

0.0 


Table 2b. Comparison of Control Forces 

Results obtained by Ref. 2 and those by PERFORM, u-. and u 4 
are taken from Ref. 2. u‘s were computed by PERFORM. 
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ment of control forces. It is a problem of suspending a 
structure within an underground cavity so as to provide opti- 
mum shock isolation. We consider the structure to be a mis- 
sile suspended in a silo in the manner shown by Fig. 10 , where 
u^ , u 2 , u^, u^, and u^ represent shock isolators to be opti- 
mized. For purposes of analysis the missile is viewed as a 
rigid body capable of three modes of in-plane motion, two 
orthogonal translations and one rotation. The silo is assumed 
to be disturbed by a shock pulse, resulting in an in-plane 
motion of the silo enclosure. This motion can be . decomposed 
into two orthogonal input motions, one vertical and one hori- 
zontal component. Consider Fig. 11 , in which the notation is 
the same as that used in Fig. 7 except that 8- = angles of 
deviation of shock isolators from the original preshock direc- 
tions. 

The equations of motion are 

4 ■ . ; : 

mz. + u c sin 8 r + 1 ' u. cos 8. = 0 ? 

1 5 p 5 ^ l l 

4 * 
mz 0 + u cos 8 C + £ u- sin 8 • = 0 

2 5 5 i=l 1 ( 17 ) 


10 

+ 

u 5 h 2 

sin 

(85 + 

0) 




. + 

U 1 L 1 

[cos 

< 8 i 

- 0) 

d 

2L 1 

sin 

(Bi - 

8)1 

- 

u 2 L 2 

[cos 

< 6 2 

- 8) 

+ d 
2L 2 

sin 

^2 

0) ] 

■ + 

u 3 L 1 

[cos 

(6 3 

+ 0) 

d 

2L 1 

sin 

(83 + 

0) ] 

-u 

4 L 2 

[cos 

<6 4 

+ e ) 

+ d : 
2L 2 

sin 

(P 4 + 

0)1 


where m is the total mass of the missile and I is the moment 
of inertia of the missile. Since 8- are functions of z, , z 2 , 

Yf, y 2 , the terms such as, 8^u^ in tfte above equations are 
nonlinear. 

Now assume, as in the previous section, that the isolator 
forces retain their horizontal or vertical directions despite 
the motions of the base and mass. Again, assume the 0 is small, 
and in addition assume ^ is small. Then the equations 
of motion are rewritten as 
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0 


mz 


z, + l u. = 
i i 

1 = l 


mz. + u c = 0 
Z D 


( 18 ) 


18 + u^L^ - + U 3 L 1 ” U 4 L 2 = 0 


The relations 


z i = x i + *1- *i * *i * »i> *i - x i + 

Z 2 = x 2 + y 2 , z 2 - i 2 + y 2 , 2 2 = x 2 + y 2 

I 

2 < 

and I = mp substituted in Eq. (16) yield, for unit mass, 

*1 + U 1 + u, 2 + u 3 +u 4 -■ -?! 

*2 *S. ", -?2 ’ (1S) 


e + U-.L - u„L + U-L - u .L_ =* 0 
1 1 2 2 3 1 4 2 


.2 . 


.2 . 


The above equations are all linear and in the desired form for 
PERFORM input format. ; 

With additional control forces u^, u 2 , u^r and u a 4 
shown in Fig. 12 and neglecting the contribution of u. , u„ , 
u,, and u. to rotational motion, the equations of motion can 
bS written as 

x. + u. + u_ + u, + u. = -y. 


y 2 + u s ■ -y 2 


( 20 ) 


0 + u .. + u» + u->+u / , = 0 
al a2 a3 a4 


Again we have equations of motion in the desired form for the 
PERFORM input format and the equations are applicable where 0 and 8. 
are not necessarily small. 
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12 Modified Silo—Missle System 


48 







C. Train Suspension Model 


The proper design of the wheel-to-car (or truck) suspen- 
sion mechanism is essential to the successful operation of a 
high speed rail passenger car. Simplified vertical and lateral 
dynamic models of a high speed railroad car have been developed 
and a trial-and-error optimization procedure has been applied 
to determine the optimum stiffness and damping properties of 
the suspension systems in Ref. 3. 

The mathematical models used were a four-degree-of-freedom 
model for vertical responses to vertical inputs and a ten- 
degree-of-freedom model for lateral responses to lateral or 
rolling inputs from the rails. The linear equations of motion 
were solved by digital computer programs for sinusoidal in- 
puts to both models and for one random input to the vertical 
model. To include nonlinearities in the lateral truck suspen- 
sion system, a real-time digital simulation computer program 
was developed and utilized. Both linear and nonlinear accel- 
eration responses to sinusoidal, transient deterministic, and 
random inputs were obtained for the lateral model by means of 
this simulation program. 

In this work we formulate the application of the computer 
system PERFORM in evaluating the limiting performance of the 
lateral train model, subject to transient inputs. 

1 . Description of the System 

The lateral train model used in Ref. 3 is a ten degree- 
of-freedom system as shown in Fig. 13. The motion of this 
model is described by the following generalized coordinates: 

- car lateral bending 

- car lateral rigid-body translation 

- car rigid body roll 

- car torsion 

- car rigid-body yaw 

- transformer lateral translation 

- rear traction motor lateral translation 

- forward traction motor lateral translation 


b 

(t) 

*c 

(t) 

6 co 

(t) 

0 cl 

(t) 

n 

(t) 

y *g 

(t) 

y Ar 

(t) 

y Af 

(t) 
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0 Ar (t) “ rear traction motor roll 


®Af(t) - forward traction motor roll 


For a limiting performance study we replace the spring 
and damping forces in Fig. 13 by the control or isolator 
forces as shown in Figs. 14 and 15. Fig. 14 represents a 
general case where the entire spring and damping forces of 
the system are replaced by eleven isolator forces. A 
special case of only 3 isolator forces is shown in Fig. 15. 
The eleven isolator forces considered are: 

u - rear lateral bolster force 
cr 

u r. - forward lateral bolster force 
cf 

u - rear vertical bolster force 

SIT -i 

u - - forward vertical bolster force 
sf ! 


u - transformer lateral bolster force 

g 


u 


Hr 


rear traction motor rolling torque 


u Af ~ forward traction motor rolling torque 
u Ar - rear vertical equalizer force 
u A £ - forward vertical equalizer force 

. . i 

u ar - rear traction motor rolling torque 
u a £ - forward traction motor rolling torque 


2 . Operating Conditions and Objectives 

The railroad car system is subjected to transient lateral 
and cross-level inputs which simulate the effects of a lateral 
deviation of the track from a straight course and the bank on 
a curve, respectively. 

The objective is to find the isolator forces that minimize 
a maximum lateral acceleration of the car. Furthermore, the 
constraints are on the lateral deflection of the bolster springs 
Thus we wish to find the unknown isolator forces such that 
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Fig. 13 Mathematical Lateral Model of Railroad 
(a) Top and Mid Section Views 


5] 


C/l£. 


B 


I 

8 




— ycC-fc) + v(T, t) 




a 'd 




Ur 


■4f.* 

- 


Car and Truck Suspension 




Fig. 13 Concluded 

(b) Truck Suspension (Section B-B, and C-C) 
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Fig. 14 Schematic Vlev of Lateral Rail Car Model 
With Eleven Isolator Forces 
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max | A . (x) 


is minimized, and 


YL<AY ( t ) < YU 
— r — 


YL<AY f (t ) <YU 


are satisfied. Here 


A cJl (x) = B (t)W(x) + ? c - h 4 (0 co +.0 cl cos ~) - (| - x)ri 
= + q 2 - h 4 q 3 - (h 4 COS ^)g 4 - - x >q 5 


2r Ar 


~ h 3 [Q co (t) + Q cl {t) 005 L“ ] j 
Ay = same as Ay with W(d) replaced by W(L-d), 

II ! IT 

,L - . . ,L .j % ird i tt (L-dy 

(2 - d > by- (3 " d) ' — • ' l ~ 

and the subscript r by f 


3* Modified Equations of Motion for PERFORM 

The equations of motion for the original system shown in 
Pig. 13 were derived in Ref. 3 and expressed in the form 
• • • 

Mq + C q + Kq = Ff fc (21) 
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In order to use PERFORM the equations of motion must appear 
as 

•• • 

Mq+Cq +Kq+Uu=Ff k (22) 

where the u’s are the control forces that replace the spring) 
and damping forces. 


Formulation of equations of motion in the form of Eq. 

(22) is not unique and depends on the choice of the coordi- 
nate system as well as the isolator forces to be used. This 
will be shown in the following illustrations with a two-degree- 
of-freedom system with its kinematic relations as indicated in 
Fig. 16a. The equations of motion in the coordinates and 
are written as 


m l^l + -^1 (Zf-'f) + k 2 ( z x — z 2 ) = 0 
m 2 z 2 + k 2 (z 2 - Zl ) = 0 ■ 

Using the coordinates and x 2 the equations of 

m l^l + k l x l ” k 2 X 2 
m 2 .« 2 +« x ) + k 2 x 2 = -m 2 f 


(23) 

motion become 


(24) 


Now we write the equations of motion in the form of Eqs. (23) 
and (24) for the following six different cases of coordinate 
systems and isolator forces . 


Case 1 (Fig. 16b) 

m l^l + U 1 + u 2 5=0 

(25) 

“ 2*2 ‘ u 2 = 0 


where 

U 1 = k l (z i -f ) 


u 2 = k 2 (z l' z 2> 


Case 2 (Fig. 16c) 


m 1 x 1 + u x + u 2 =-m x f 
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m 2 (x 2 + x^ ) -u 2 


(26) 


“ -m 2 f 


where 


u i = k i x i 


u 2 = ' k 2 x 2 


Case 3 (Fig. 16d) 

m l^l + U 1 + u 2 = 

m 2 S 2 -u 2 - 0 
Case 4 (Fig. 16e) 

m l^l + k l Z 1 + u 2 = k l f 

m 2 2 2 - u 2 = 0 
where 


(27) 


(28) 


U 2 k 2 (z l z 2 ) 


Case 5 (Fig. 16f) 


m l^l + k l X 1 + u 2 = 
m 2 (x 1 + x 2 ) - u 2 = -m 2 f 

where 


u. 


= -k. 


Case 6 (Fig. 16g) 


m l x l + k l x l + u 2 


-m^f 


m 2 z 2 " u 2 = 0 


(29) 


(30) 


In Case 1, where the coordinates and z 2 represent absolute 
displacements of masses and m 2 ,the equations of motion 
(Eq. (25)) do not have any forcing functions. This would not 
be an acceptable form unless only initial conditions of these 
coordinates are specified. Furthermore, if for Case 1 we 
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define 


u^ = k^ and 

u 2 = k 2 ^ Z 1 ” z 2 ^ ' 
the equations of motion become 

m i h + u i + u 2 = k i f 

(31) 

In 2 S 2 ' U 2 = 0 


Eq. (31) includes a forcing function and may appear to be as 
acceptable as Case 4. However, we have here the problem of 
a term with the system spring characteristics with k-^ not 
being completely replaced by u^ or u 2 . Cases 2 through 6 
are all acceptable and the choice of a particular coordinate 
system or isolator forces depends on the initial conditions 
and the nature of the problem at hand. 


a. The Case of Eleven Isolator Forces (Fig. 14) . 


For this case the equations of motion in the form of Eq. 
(22) will have the following vectors and coefficient matrices: 


. Displacement Vector 
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where 


z 

r 

II 

> 

h 

-Y 

r 

z f 

ii 

> 

Hi 

- Y f ■ 

^r 

II 

CD 

> 

H 

5 r 

R 

*f 

- e Af 

5 f 

R 


The coordinate system used here is similar to that of Case 3 
discussed above. . 


Forcing Function Vector 

r**~ 

f ll 


f l * 


12 


13 


14 


s . 


where 


f,, = M. Y 
11 Ar r 


f = M Y 
12 Af f 


f . 6 

13 Ar r 


f 14 I Af 6 f 


and 

Y of 


2 ttV ( t — t ) 
o 


f f _ 2 [1 “ C ° S X + (A£) t 3 ' 


2 ttV ( t~ 1 2 ) 


Y r “ 2 11 Cos X + (AA) t 3 ' 


Y f = 


2 Y Of (TTV) 2 

[X + ( A £ ) t ] 


Y r= 2Y or UV) 


IX + (AA) t ] 


Cos 


Cos 


2W (t-t Q ) 
X + ( A £ ) , 


2 ttV (t-t 2 ) 

X + TEXT 
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[1 - Cos 


Cos 


6 = 
r 


5 f - 


or 


of 


2ttV (t-t ) 
o 

A + (A£) . 


1 , 


[1 - Cos 


2irV (t-t 2 ) 

r~rrnT7 


6 = 
r 


2 6 (ttV) 
or 


[A + CA£) t J 


26 


6 f =' 


of 


(ttV) 


[A + (AA)^] 


Cos 


2ttV (t-t ) 
o 

A + (A£) 7 


2ttV (t-t 2 ) 
A + (Ail) 7 


The Y's are lateral inputs and the 6 ' s are the cross-level 
inputs. The second derivatives of Y's and 6's with respect 
to time are expressed as Y's and 6's. The quantities V, A, 
and, : (Ail) t represent the train speed, the wave, length of the 
input disturbance, and the wheel base of a truck suspension, 
respectively. In the above equations is obtained from the 

A+ ( A £ ) . ' . 


relations 


t, - t 
•1 o 


fc 2 ' fc l = 


V 

L - 2d - (Ail) - A 
V 


where t. - t 
1 o 


= time interval for forward truck to roll over 
the disturbance. 


t 2 ~ t, = time interval between when the front truck 
leaves the disturbance and the rear truck 
encounters it. 

. Isolator Force Vector. 


u = [u cr = u x , u cf = u 2 , u sr = u 3 , 
u s f = u 4 , u g = u 5 , u ir = u 6 , 
uilf = u 7 , UAr “ u 8 , UAf = u 9 , 
u ar = u ]_q ' u af = u n^ • 



where the isolator forces are defined in terms of the confi- 
guration of Fig. 13. They are 


ird 


u = u. = c 0 [W(d)£> + y .+ h. 6 + h_ cos T „ . 

cr 1 x-cr J c 3 co 3 L cl 2 


6 


- - d)n 


-d r + 4 r > + h 2r (* r + X-)] + k lor [W(d)..£> ;+ y c + h 3 0 


CO 


+ cos tt © c i - (t ~ d )n - ( z r + + h ov (<Pr + I 5 -)] 


2 r 


u f = u Same as u with the subscript r replaced by f, 
C W(d) by W^E-d) , 

,L ,, , ,L , . . nd , ir(L-d) 

-(^ -d) by -d) , and by — ^ 

• * 

U =U =C +x COS ~ 0 . - Y (<P + ■—) ] 

sr 3 s 2 co 4 L cl 4 Y r R 


ird 


+ k [±- 8 + t cos ”T~ 6 , - t (<P + 

s 2 co 4 L cl 4 Y r 


■)] 


u f = u = Same as u with d replaced by (L-d) and the sub- 
'~' L script r Sy f. 

= U c = c n „ [W {^) b + y„ + h„ 0 

,L 


5 ~£g l " '2' " ‘ y c + n T 0 co Y g J 


+ k £g [W ¥ b + y c + h T 0 

6 _ 


- Y ] 
co g 


u_ = u c = c„ [z + h. (d> + - 7 T-) ] + k. [z + h-. (d> + - 5 -)] 

6 £r x lr Y r R lx x lr Y r R 


ir 


u^ f = u^ = Same as u^ r with subscript r replaced by f. 

& 6 

u, = u Q = c. (d> + + k,. (d> + 

Ar 8 Ar Y r R Ar Y r R 

u^ = Ug = Same as u^ r with subscript r replaced by f. 

6 6 

u ar - u 10 ’ c ar + IT’ + k ar <*r + TT 1 

u a ^ = u^^ = Same as u ar with subscript r replaced by f. 

. Damping Matrix 
C = [0] 
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• Matrix 

- 18 foj 

exc ept 


K. 


II ® (Bl) 


■-Z& 


K ^"¥, 
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The Case of Three Isolator Forces (Fig. 15) 


The vectors and coefficient matrices of the equations 
of motion in the form of Eq. (22) are defined as follows: 

. Displacement Vector 



Here we are using a coordinate system similar to that of 
Case 4. 

. Forcing Function Vector 



where 
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f ll 

= k ir 

Y 

r 

+ 

c £r 

Y 

r 

f , 

= k 

Y 

+ 

c „ „ 

• 

Y„ 

12 

Af 

f 


If 

f 

f , „ 

= k. 

6 

+ 

c , 

• 

5 

13 

Ar 

r 


Ar 

r 

£ 14 

4-) 

< 

X 

II 

S £ 

+ 

C Af 



The rail inputs Y's and 6's are the same as in the first 
case. Y's and 6's are the time derivatives 


Y = 
r 


Y irV 
or 

+ (AH) 


2ttV (t-t 2 ) 

x + (a rr ' 


Y f TTV 

Y = — °^ 

f X + (A £) 


2irV (t-t ) 
o 

X + (A£) , 


• • 

The expressions for 6 and are of the same form as for the 

Y and Y,. with 6 - anS 6 replacing Y _ and Y 

r f of or ^ ^ of or. 


f of or 

. Control Force Vector 


U cr = U 1 

5 = U cf = U 2 

u = U-, 

g 3 

where the isolator forces are defined in terms of the 
configuration of Fig. 13. These are 

U 1 = c Hcr [W|d)6 + + h 3 9 co + h 3 COS TT 9 cl "'I - d) * 

- ?Ar + h 2r e Ar! 

+ k lcr [W(d)b + y c + h 3 e co + h 3 cos ^ 0 ol - (| - d)n 


_ y Ar + h 2r 0 Ar"' 

u 2 = Same as u^ with the subscript r replaced by f, and W(d) 
by W(L-d),-(| -d) by <| - d) , and ^ by 
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U 3 " C W [W( I> i + + h T S co ‘ V + k 2g [M( ^ b + h n 

+ Yc - yg] 


. Mass Matrix 

M = same as in the previous case of 11 isolator 
forces 

. ..Stiffness Matrix, K 


K 11 

— 

(El ) L 

K 12 

= 

... = K 

-L / 

K 21 

= 

• ■ • — ' 

^ / 

K 31 

= 

K 32 = 0 



k B 2 

K 33 

— 

s 

2 



k B 2 

K 34 

— 

4 [cos 

K 35 

= 

... = K 38 


= 0 
= 0 


•rrd 


, tt (L - d) , T , 

L + COS L ] = K 43 


-k s B2_ 


K 39- — — - K 93 = K 3 ,i 0 = Kio,3= K 49 = K 4 ,l 0 = K 10 , 4 


K = 

■ 44 4 

r 2 

- [cos 

• 

• 

• 

II 

m 

II 

•fe. 

00 

II 

• 

• 

• 

II 

rH 

in 

* 

K 5,10 

• 

• 

• 

II 

iH 

VO 

ii 

w 

cn 

*•» 

M 

O 

M 

II 

• 

• 

• 

ii 

« 

OV 

II 

* 

-J 

II 

• 

• 

• 

= k • 
lr 


2 L 


= 0 
= 0 


co 


= k 94 
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0 


K 78 K 87 


K 88 = k £f' 


K 89 = K 98 =0 


k 8,10 = K 10,8 * k A fh lf 


K 


K 


91 

K 92 

K 95 K 96 

9 7 

l) 

>3 

-4 

11 

k fr h lr 


k B 2 

K r2 


s 

+ Ar f 

99 

4 

+ 4 + 


" 9,10 


‘ar 


lr R £r 


'44 $ 


ks 2 K a£ R 2 2 

K 10,10 = 4 + 4 + k af + k Af h lf 

Damping Matrix, C 

Identical in form to K except = 0, 

c s B 2 , 0 rrd 

[cos* + cos 
L 


2 fr(L - d) 


For other elements k , k. , k a are replaced by c 


C A' c a respectively. 


, Control Force Matrix, U 
U il * W < d > 

U 12 
U 13 
21 
? 31 
J 33- 


U 
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~d 


U 


W (L - d) 

U 53 

~ °6l 

’ U 62 

= 0 

»<!> 

■ ; u 6 3 

U 71 

= -1 


U 22 = U 23 ~ 1 

U 72 

“ U 73 

* °81 

* 0 

U 32 = h 3 

U 82 

= -1 


! 

h T 

- U 83 ■ 

= 0 
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u 


41 


Cos 


ird 

L 




U 


42 


h 3 Cos 



U 


92 


= U 


93 


~ U 10,l ~ U 10 , 3 




U 


10,2 




. Forcing 

Function 

Matrix 



(1) 

(2) 

(3) 

(4) 

(1) 

0 

0 

0 

0 

(2) 

0 

0 

0 

0 

(3) 

0 

0 

0 

0 

(4) 

0 

0 

0 

0 

(5) 

0 

0 

0 

0 

(6) 

0 

0 

0 

0 

(7) 

1 

0 

0 

0 

(8) 

0 

1 

0 

0 

(9) 

h, 

lr 

0 

R 

4 

0 

(10) 

0 

h lf 

0 

R 

4 _ 


4. Notation 


A c« <x) 


lateral acceleration of car 


B 

b (t) 


lateral distance between bolster springs 
(see Fig. 14) 

-generalized coordinate for lateral car 
bending 

(h-1) ^element of damping matrix 


d 


distance from end of car to centerline of 
trucks (see Fig. 13) 
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Young's modulus 

0) 

frequency of qth degree of freedom, Hz 


modulus of rigidity . 

vertical distance from traction-motor 
center of gravity to line of action of 
lateral equalizer spring (see Fig . 13 (b) ) 

vertical distance from traction-motor 
center of gravity to line of action of 
lateral bolster spring (see Fig. 13(b)) 

vertical distance between car elastic axis 
and line of action of lateral bolster 
spring (see Fig. 13(b)) 

vertical distance of car section center 
of gravity from car elastic axis, positive 
for center of gravity above elastic axis 
(see Fig. 13 (b) ) 

vertical distance between car elastic 
axis and line of action of lateral trans- 
former spring (see Fig. 13(a)) 

flexural moment of inertia of car cross- 
section 

mass moment of inertia of traction motor 
in roll about its center of gravity 

mass moment of inertia of car in roll 
about car elastic axis , per unit length 

torsional constant of car cross-section 

spring constant 

(h-i)th element of stiffness matrix 

length of car (see Figs. 13(a)) 

(h-i) t * 1 element of mass matrix 


mass of car per unit length 
traction-motor mass 



M 


C 



R 

S (t) 


v 

V (x,t) 


W(x) 


X 


Y (t) 

YL . 
YU 

V t} 

y c (t) 
6 (t) 


(A*) t 


car mass, mL 
transformer mass 


viscous damping coefficient for 
degree of freedom, where q = A, cl, cT, 
etc., also see subscripts 

i 

lateral distance between equalizer springs 
(see Fig. 13 (b) ) 

vertical displacement forcing function at 
rails (see Fig.. 13 (b) ) 

car speed 

lateral bending deformation of car (see 
Fig. 13 (a) ) 

car bending mode-shape value at x 

length along car, measured from rear of 
car 

lateral displacement forcing function at 
rails (see Fig. 13(b)) 

lower limit 

upper limit 

generalized coordinate of lateral trans- 
lation of traction motor 

generalized coordinate of lateral rigid- 
body translation of car (see Fig. 13(b)) 

generalized coordinate of lateral trans- 
former translation 

increment of vertical displacement at rails 
associated with cross-level (or rolling) 
displacement input (see Fig. 13(b)) 

center- to-center distance between axles 
of a truck suspension system 


Ay (t) 


displacement relation in lateral bolster 
springs for nonlinear spring behavior 
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n (t ) 


generalized coordinate of car in rigid- 
body yaw 


e A (t) 


generalized coordinate of traction-motor 
rolling degree of freedom 


6 C (*,t) 


rigid-body roll and .torsion of car 


6 


co 


(t) 



(t) 


Subscripts 

A 


c 


generalized coordinate for rigid-car roll 

s t 

generalized coordinate for 1 torsion 
mode of car 


identifies properties associated with the 
traction-motor-equalizer system 

associated with car 


f 


identifies forward trucks 


g 

h, i 

a 


He 


associated with transformer 

integers identifying matrix elements 

identifies properties of traction-motor- 
lateral-equalizer system 

identifies properties of lateral bolster 
spring-damper system 


o 


identifies amplitude of forcing function 


r identifies rear trucks 

s identifies properties of vertical bolster 

spring-damper system 

Dots over quantities denote differentiation with respect to 
time . 
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D . Aircraft Ride Control System 


To evaluate the candidate designs for the ride control 
system in a STOL airplane, a limiting performance problem 
was formulated. Specifically, it is desired to minimize 
the maximum vertical center-of -gravity (c.g.) acceleration 
subject to constraints on flap and elevator deflections. These 
two deflections are used as controls. 


The airplane is assumed to behave as a rigid body. The 
nomenclature is: 

a = wing angle of attack 

0 = rigid-body pitch of airplane 

v = perturbation in forward velocity 

a (t) = angle of attack due to gust 

oig(t-T) = delayed gust angle of attack (at tail) 

T = time delay between impingement of gust on 
wing and tail 

u^ = elevator angular deflection, a control function 

u 2 = flap angular deflection, a control function 

a 0 = gust amplitude 
a) = gust frequency 

Assuming quasi-steady aerodynamics, the equations of 
motion for the airplace in first order PERFORM format are 

s — A s + Bu + D f^ (32) 

where 


- — 



S 1 


8 

S 2 


• 

8 

S 3 


a 

S 4 


V 
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u 


(34) 



o g (t) 
a (t - T) 

. g 

Consider the specific case: 



A = 


0 

0 

0 

-4.77992 


1.0 

-0.05662 

0.97710 

0 


0 

-0.58024 

-1.75148 

0.96138 


(35) 


0 

-0.0000262 
- 0 . 001 - .. 

0 


B 


0 

-0.64961 

-0.15183 

0 


0 

-0.18396 

-0.64855 

0 


D 


0 

-0.13587 

-1.65057 

0 


0 

-0.10674 

-0.14425 

0 


Select VqIs^ - SO as the objective function, with v 0 = 

308 ft/sec. The constraints are imposed on the magnitude of 
the control. Numerical solutions are obtained for two cases. 
In Case 1 only u^ is active and U2 = 0. In Case 2, both u 
and Uo are active. The disturbance in both instances is a 1 
step input, i.e., a g (t) = 0.649H(t) where H(t) is a unit step 
function. The time delay is T = .08207 sec. The total time 
interval used was 0.95 sec with a time increment of 0.05 sec. 
The trade-off curves for these two cases are shown in Fig. 17 . 
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E. Shock Absorbers for Freight Cars 

Consider the limiting performance of the shock-absorber 
system (or cushion) of railroad vehicles used in protecting 
passengers or cargos under crash conditions. Typical analyses 
of freight car lading protection problems are described 
in Ref s . 4 and 5 . 

The dynamic system model of Fig. 18 is taken from Ref. 5, 
with the following nomenclature: 

S Q = lading contact area 

c = damping coefficient, due to friction in lading 

k = elastic constant, representing resilience in lading 

’ • ' } 

M c = mass of struck car body 

■ , t 

, • ‘ *• t £ “*'■ 

M = mass of lading s 

±J 

- i z.~ 

= mass of striking car; . t j 

u = cushion force 
d = y-L “ y 2 = cushion travel 

V = impact velocity 

.j... • V'* * -*.v % 

The equations of motion are given by . 

My. + u = 0 

S'* 1 

M c y 2 - k(y 3 - y 2 ) - c (y 3 - y 2 > - u = 0 (36) 

Va + k( y 3 ' y 2> • c<y 3 ■ y 2> = 0 
We will treat the impact conditions 
y x (0) = y 2 (0) = y 2 (0) = 0 

y 2 (0) = V 

where t = 0 is the time at which impact occurs. 
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Fig 18. Model of striking car, struck car and 
la d ing. 



For specified impact velocities, the problem is to mini- 
mize the peak force transmitted to the lading while the cushion 
travel distance is bounded. Thus, the u is sought that 
minimizes 


- max | k (y 3 - y 2 ) i (37) 

with the restriction that 

i* 

0 1 d <_ A (38) 

where A is prescribed. 

The problem has been put in PERFORM format for the fol- 
lowing physical and material constraints; 

c = 23,000 lb sec/ft 

k = 341,667 lb/ft > • 

M = 1630 lb sec 2 /f t \ : . 

= 1590 lb sec 2 /ft 

M = 5280 lb sec 2 /ft 
s 

S = 33.7 ft 2 • ' ' 

o 

i 

V = 10 mph , 

A PERFORM generated trade-off diagram between maximum 
transmitted force and cushion travel distance is shown in 
Fig. 19. 
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Ltted to lading {lb/ft ) 

|k(y 3 -y 2 ) |/s 0 



0> : . 10 20 30 : 40 

cushion travel distance, A (inches) 


Fig 19, Trade-off diagram for rail vehicle cushioning system 
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P. A Launch Booster Control Problem 


We consider here a booster control problem known as the 
constrained bending moment minimax drift problem. This is the 
problem of minimizing the lateral drift from a reference 
trajectory along a flight path where the magnitude of the bend- 
ing moment cannot exceed a critical value (Refs. 6 , 7 ). This 
problem will be placed in the PERFORM format in the following 
subsections . 

1. Equations of Motion 

The rigid body equations of a typical launch booster are 
given in the following form (Ref. 6 : 

(39) 

(40) 

(41) 

where <J>, a, 3, z represent angular deviation from reference, 
angle of attack, engine gimble angle, and lateral positional 
deviation from reference, respectively. The wind contribu- 
tion to the angle of attack is denoted by a y and the vehicle 
speed by V. We assume several transient wind disturbance 
profiles which might be encountered by the launch vehicle 
during flight are known. 

By substituting Eq. (41) into Eqs . (39) and (40), we 

have 

♦ = 'V + ' -T' Z - C l% ' C 2 6 

K 

* - < K 1 + K 2»* - IT i * Va + K 3 6 

Let 


• 



<p 


rs i 

• 



z 


s 2 

<J> 


S 3 

z 


S 4 

— . 


_ 


<j> + + C 2 3 = 0 

Z = + K 2 ct + k 3 3 

a = <p - ^ + a 
Y V a) 
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and u = 3 / and f . = a 

k 03 


Then the rigid body equations of motion for the booster take 
the PERFORM compatible form 


s=As + Bu + Df^ 

where 



— 

C. 


— 


0 

1 

V 

- C 1 

0 


0 

- K 2 

V 

K l +K 2 

0 

A = 






1 

0 

0 

0 


_ 0 

1 

0 

0 


_ - C 2~ 



■- c r 

B = 

K 3 

and 

D = 

K 2 


0 



0 


_ 0 _ 



_ 0 


2 . Objective Function and Constraints 

The bending moment, given by a linear combination of a 
and 3 is constrained as follows: 


I M' a + JA' a 3 I < M 
1 : a 3 ' — 


max 


(42) 


where M^ and M^g are bending moment coefficients and Mmax i s 
the maximum allowable bending moment. We choose as an objective 
the minimization of the maximum | Z (t) | of the booster subject 
to the constraint given by Eq. (42). 

In the PERFORM format the objective function is written 
as 


PXl s + PX2u + PX3 f 

where 

PX1 =[0100] 
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PX2 = [ 0 ] 

PX3 = [ 0 ] 

J 

In the PERFORM input format the constraint given by Eq. 
( 42 ) can be expressed as 

YL < Yls + Y2u + Y3 f < YU 


where 


YL 



YU 

Y1 


' M max ‘ 

■ K 

to — rr 


M a 0] 


Y2 =, [M" 0 ] 


Y3 = ! [M" ] 
— a 



The specification matrix is given by' 

MSP =[003 

) 

i 

3 . Output 

The rate of ^change of the lateral positional deviation 
from reference, S(t), can be tabulated versus time. Also, 
a trade-off relation between the minmax Z(t). and the maximum 
bending moment can be obtained by varying the bounds of the 
bending moment. 
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G. 


A Pressurized-Water Nuclear Reactor Power Plant 


As a potential application of the computer program to 
a large system, we consider the control of a pressurized- 
water nuclear power plant. This system was first used as an 
example in Ref. [ 8] . Fig. 20 shows the schematic diagram of 
the pressurized water reactor. Heat is produced in the reac- 
tor by nuclear fission, with the rate controlled by insertion 
and withdrawal of control rods. This heat is carried away by 
the primary coolant which then flows to a heat exchanger where 
it gives up its heat to make steam. The steam goes to a tur- 
bine that converts the energy to shaft work to produce elec- 
tric power. Pumps cause the primary coolant, to circulate 
through the primary loop. 

A pressurizer rides above the primary loop with saturated 
water in the lower part and saturated steam in the upper. 
Having saturated conditions in the pressurizer assures that 
the water in the primary loop is subcooled. The steam in the 
upper part of the pressurizer acts as a "cushion" to allow 
expansion or contraction of the water as its temperature 
varies. Spray water taken from the coldest part of the pri- 
mary loop condenses steam and lowers the pressure. An elec- 
tric heater boils the saturated water, thus lowering the 
water level and increasing pressure. 

The principal disturbance to the plant operation comes 
from sudden changes in steam demand from electrical load 
changes. The most difficult case is a load drop from maximum 
load. The reduced steam flow takes less heat from the heat 
exchanger, causing temperature and pressure to rise rapidly 
in the primary loop. 

A ten percent step drop in steam flow from full load is 
applied with the plaint initially at equilibrium. This load 
change is treated as a forcing function. 

The objective is to maximize the minimum value of the 
pressure of the pressurizer with continuous constraints on 
maximum pressure and the three controls plus terminal con- 
straints on a temperature and its rate of change. We sketch 
the formulation of this problem in the PERFORM format in the 
following subsections. The specific coefficients involved 
can be found in Ref. 8. 

1 . ' Governing Equations 

The equations describing the system are expressed as 
S=As+Bu+Df k 
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where 


s = 


u = 


f, = 


A = 


B = 


D = 


25 I 


= R 


u ' = S 

, 2 »- 


state vector 


Rod speed 
Spray 


1 U 3 = Q : 

• Heater 

? 

I f 1 

r 

L 11 J ; 


C Ai 0 

25 x 25 

L Bi 0 : 

1' , 

25 x 3 

L Dij ] 

25 x 1 


Control Forces 


Forcing Function 

‘ ■ i: •• 

Coefficient Matrix 

Coefficient Matrix 
Coefficient Matrix 


2 . Objective Function 

The objective function is expressed in the matrix form 
PX1 s + PX2 u + PX3 f, 


where 


s, u, and f, are defined above, 
k 


PX1 =) PX. .1 

L 

1 x 

25 

Coef f icient 

Matrix 

PX2 =f PX. 

13 

1 X 

3 

Coefficient 

Matrix 
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PX3 


lxl 


Coefficient Matrix 


= [ PX ii) 

3 . Constraints 
Constraints in matrix form are expressed as 
YL < Yls + Y2u + Y3f k £ YU 


where 


s, u, 
YL = 

Yl = 

Y2 = 

Y3 = 



are defined above 
6x1 

6 x 25 

6x3 

6x1 


Lower Limit 


Coefficient Matrix 


Coefficient Matrix 


Coefficient Matrix 



YU . . 6x1 


Upper- limit •• 


The specification matrix is given by 


MSP 


0 0 
0 1 
0 3 
0 0 
0 0 
0 0 


4 . Output 

The desired output vector V is given by 
V=Qls+Q2u+Q3f k 
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where 


V 


V 

V 

V 

V 

V 

V 


1 

2 

3 

4 

5 

6 


Pressure 

Temperature 

Time derivative of temperature 

Rod speed 

Spray 

Heater 


~ [ Q± jJ 6x2 

~ = [ Q± 0 6x3 

Q3 = [ Q ij] 6x1 


Coefficient Matrix 

Coefficient Matrix 
Coefficient Matrix 


Trade-off relations between the minimum_value of the pres- 
surizer pressure and the output vector V shown above can be 
obtained by varying the bounds YL and YU. 


» 
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SECTION V 


CONCLUSIONS 

A computational capability, PERFORM, for the evaluation 
of the limiting performance of transient dynamic systems has 
been developed and is described in this report. This capa- 
bility is intended to be used by a designer of dynamic systems 
to determine the feasibility of the proposed design specifi- 
cations. Thus, PERFORM is meant to be a new design tool. 

In addition to being able to use PERFORM to scrutinize 
design specifications, this capability can be used during the 
actual design process to measure the relative success of pro- 
posed designs with the characteristics of the theoretically 
best design provided by PERFORM. 

PERFORM is used in an engineering problem-oriented form. 
The user prescribes in a simple fashion the system equations 
of motion, sets of constraints, objective functions, and 
classes of possible inputs. The capability then automatically 
computes the limiting performance characteristics requested 
by the user. 
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APPENDIX I 


USERS GUIDE FOR PERFORM 


, PERFORM is a computer system that can be used to determine 
the limiting performance characteristics of a dynamic system 
described by a system of first or second order equations of motion. 
The information, including input requirements, necessary to use 
PERFORM is described in this appendix. 

A. PERFORM Utilization Procedure 

, This. section describes the steps necessary to solve a limit- 
ing performance problem with the PERFORM system. Essentially this 
is a statement of the details of the system flowchart, of Fig. 1, 
Section III. 

1. Prepare the limiting performance problem specification 
data. This input is described in detail in the , following 
sections of this appendix. In the system flowchart this 
input is labeled "PERFORM problem specification data." 

2. Run PREPROC using the data of Step 1.. There are two card 
outputs of this program; the PSTPROC Report data , and the 
linear programming problem data for the LP Solver. 

Hold the PSTPROC report data output for use in Step 6. 

'If several sets of constraints are being run for trade- 
off relations, one LP Solver data deck will be produced 
for each different constraint. . 

Deck preparation for compiling and executing' the PREPROC 
program is shown in Fig. 1-1. The details of the leading 
control cards depend on the operating system^of . the 
. computer and any user desired system optipn. 

3. Transmit the linear programming program data to the 
external computer, i.e. the computer with either MPS/360 
or a MPS/360 compatible linear programming software, for 
solution of the linear programming problem. 


If recurring use of PERFORM were anticipated, excessive recompila- 
tion can be avoided by storing the program, in object form, on the 
system library. The indicated deck structure would then have to 
be appropriately modified. 
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4. Run MPS/360, or compatible software, to obtain the 
problem solution. For each constraint deck produced 
one execution of MPS is required; however, only one job 
and transmission is required. 

The deck contents for an MPS/360 run is shown in Fig. 1-2. 
.'This model, in general, would be correct for any MPS/360 
installation using OS/360. Words in upper case letters 
are invariant among jobs and installations. Words in 
lower case letters depend on the job and installation. 

Those dependent on the job are user chosen to identify 
the problem. Those dependent on the installation should 
be obtained from system personnel. 

Because. of possible idiosyncrasies of the external computer, 
. consultation* with system personnel before attempted use is 
" imperative. ‘ 

5. ^ Receive from the external computer the linear programming 

problem solution. This can be either in card form via 
. courier or by telecommunications if facilities are 
available. 

6. Run PSTPROC using the MPS/360 linear programming solution 

data received in Step 5 and the PSTPROC report data of 
Step 1 . . 

Fig. 1-3 shows the deck structure for this run. The de- 
tails of the leading control cards depend on the operating 
system being used and any user desired system options.* 

The different linear programming solution decks must be 
physically concentrated behind the PSTPROC data. The 
importance of careful card handling cannot be understated 
since the existing PERFORM system has no facilities for 
error detection. 

The limiting performance problem solution is output in 
this step. 



As; indicated in the discussion of Step 2, it may be advisable to 
place the PSTPROC program on the system library. 
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Job and other Control Cards 

T ' ‘ ‘ ’ 

PREPROC program 


Data Indicator (CDC 6000.: 7-8-9. card) 

■j ' ' ' 

Problem Specification Data 


End of Job Indicator (CDC 6000: 6-7-8-! 


(Lines Represent Decks of Cards) 


Figure 1-1 PREPROC Deck Structure. 
Compile and Execute 


card) 
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account data 


// Job-name 

// 


JOB 

EX 


MPS- procedure -name 


// Compile-step- name. SYSIN 


t 


/* 


// Go-'s t epriame . S Y s IN 


DD 


MPS'/ 3 60 CONTROL PROGRAM 


. • .:/• i *' V 


DD 




MPS/360 LINEAR PROGRAMMING DATA 


/ i 


(Lines Represent Decks of Card) 


Figure 1-2 MPS/360 

* t 


Execution Deck 
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Job and other Control Cards 

* 


PSTPROC program 


Data Indicator (CDC 6000: 7-8-9 card) 

A* 


PSTPROC Report Deck 

t 5 

MPS/360 SOLUTION Deck 
| (REPEATED) 


End of job indicator (CDC 6000: 6-7-8-9 card) 

i 


(Lines Represent Decks of Card) 


Figure 1-3 PSTPROC Deck Structure 
Compile and Execute 
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B. Input for PERFORM 


1 . System Described by Second Order Equations 

PERFORM treats systems described by the second 
order equations of motion 


n n n nu nf 

l + l Ci-A + l K i1 q. + l U^u = l F f (X-l) 

j=l J j=l 3 j=l 13 3 j=l 13 3 j=l 13 K3 


i = 1 to n 

in matrix form these become 


Mq + Cq + Kq + Uu = Ff 


( 1 - 2 ) 


where 


q = 


displacement vector 


n I 


n x 1 


*n 


n x 1 


velocity vector 
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r *• 


q = M3' 


l n 


n x 1 


acceleration vector 


u = 


u. 


u. 


u 


nu 


nu x 1 


control force vector 


V 


kl 


k2 


"knf 


nf x 1 


forcing function vector. The 
subscript k identifies the set 
of forcing functions, e.g. f 21 

is the first element of the 
second set of forcing functions, 
This allows the system to 
encounter alternate sets of 
forcing functions which might 
occur with equal probability. 


M = 


11 

M 

12 

M ln 



21 

M 22 

*'* M 2n 

- fMijJn x n 

mass matrix 

nl 

M n2 

... M 

nn 
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c 

K 

U 


[C. .] 
l] n x 

[K. .] 

1 } n x 

[U. .] 

30 n x 


n 

n 

nu 


damping matrix 

spring matrix 

matrix of constants 
associated with control force 


F 


[ F . ■ J r 

13 n x nl 


matrix of constants 
associated with external 
disturbances 


The user should place his equations of motion in the form 
of Eq. (1-2) . The non-zero elements of the matrices M, C, K, U, 

F are then entered as inputs in the form explained subsequently. 
This is accomplished by identifying the matrix, e.g., M MATRIX, 
and then specifying an 1 element and its value, e.g., i, j and M... 
Elements not entered are assumed to be zero. ^ 


Example 1 

The equations of motion for a system shown in Fig. 1-4 can be 
written as 


m^x + k(x-y) + u = n^g + h k (t) 


m 2 y + cy + k(y-x) = m 2 g 
Put in matrix form 


— 

— 


— — 


— 



— — 

m. 

0 


X 


0 

0 


• 

X 

1 




+ 





0 

m 2_ 


y 


0 

c 


• 

y 




(1-3) 


Now define 


x = q x y = q 2 


q = 


X 

Y 


m. 

1 


1 



m 2 

0 


— 

I L 


g 

h k (t) 


(1-4) 
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Note that the k subscript of f designates one set of forcing 
functions. In our case f^ = g = gravitational acceleration 

constant while f ^ = h^(t) = some function of time, e.g. those 

in Fig. 1-5. If f ^ 2 ’ ^ 22 ' an( ^ f 32 of F; *-9* I- 5 were all to be 

included in the input, then these disturbances would be treated 
as though all are potential forcing functions that may occur with 
equal probability. Using Eqs. (1-5), Eq. (1-4) becomes 


101 



2. 


•« • 

M q + C q + K q + U u = J 


■t. k ; 


( 1 - 6 ) 


System Described by First Order Equations 

The first order equations of motion accepted by 
the computer system are expressed as 


s = As + Bu + Df^ 


( 1 - 7 ) 


where 


s = 


m 


m 


m x 1 


m x 1 


state vector ! c-'V *• 


time derivative of state vector 


A 

B 

D 

m 


[A. .] 

lj m x m 


coefficient matrix 
coefficient matrix 
coefficient matrix 
number of equations m the system of equations 


[B. .] 

13 m x nu 

[D. .] , 

lj m x nf 


The user should input non-zero elements of matrices A, B, D 
in the form explained later. This is accomplished by identifying 
the matrices; e.g., A MATRIX, and then specifying an element and 
its value, e.g., i, j and A(i,j). Elements not entered are assumed 
to be zero. 
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Example 2 

Equations of motion for the system shown in Fig. I 
be written as 

dh. h.-h- 

c rai =— ir 1 + e(t > 


dh_ h -h, h- 

c 2~at = sf + or 

-i i i 

~dt = R^jC^ h l + R^ h 2 + C 1 e(t) 


dh 


« ' rJc 2 h i - ( r 2 c 2 + R 2 C 2 )h 2 + C 2 E k (t) 


The system of Eqs. (1-8) is equivalent to 
s = As + Bu + Df k 


where 


s = 


u = [e (t) ] 

V iE k (tn 


A = 


R 1 C 1 


R 1 C 1 


R i C l 


-(.A_ + -I_) 
R 1 C 2 R 2 C 2 


C. 


B = 


-6 can 


(1-8) 


(1-9) 
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h^, ^2 = level in tank 1,2 respectively 

C^, C 2 *» cross sectional area of tank 1, 2 respectively 

R^, R 2 = orifice resistances 

Flow rate through orifice «* — (assumption) 

K 


Fig. 1-6 Two-tank liquid-level system 




0 


D 



3 . Objective Function 


Regardless of the form (first or second order) used to 
describe the equations of motion, the format for the objec- 
tive function is the same. In the case of the second order 
equations, the user should establish a state variable vector 
s as F - 1 



( 1 - 10 ) 


q 


The user can choose any linear combination of state 
variables, derivatives of state variables, or control forces 
as an objective function. In the case of the system des- 
cribed by second order equations, these become linear com- 
binations of accelerations, velocities, displacements and 
control forces. PERFORM finds the vector u that minimizes 

the maximum (ISP = 0) (or maximizes the minimum (ISP = 1)) 
time values of the functions used as objective functions 
while satisfying the constraints. The objective function 
is formed as 


PXl s + PX2u + PX3f j. 


( 1 - 11 ) 


If more than one row of the matrices in Eg. (1-11) contains 
non- zero elements, then the peak values in time of the 
vectors resulting from the meaningful rows are to be com- 
pared. PERFORM minimizes (maximizes) the maximum (minimum) 
of the peak values. 


Example 3 

To illustrate the formation of an objective function, consider 
the 2-degree-of -freedom (DOF) system in Example 1. The equations 
of motion are given as Eqs. (1-3). Suppose we have only one set of 
forcing functions. We assume the objective function is to minimize 
the maximum of the peak values of jx]; |y|; - $\ . From Eqs. (1-3) 

y and x can be expressed as 
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-k k . + h k (t) 

x = — x + — y - u + g + 

m 1 m l m l 


y = -py + ^— x-^-y + g 
m 2 m 2 m 2 


Then the objective function involves 


k k h k ^ 

— x + — y-u + g + 

m^ m l m l 


- — y + — x - — y + g 
m 2 m 2 m 2 


( 1 - 12 ) 


I * ‘ * I 

PERFORM will minimize the maximum of peak values of these three 
relations. In terms of the form required for PERFORM, we write 



u = [u] 

g 

V 

\ (t) 


so that the objective function becomes 


0 

0 

-k 

m l 

k 

m l 

s + 

-1 

u 

+ 

1 

1 

m l 

0 

- c 
m 2 

k 

m 2 

-k 

m 2 


0 



1 

0 

1 

-1 

0 

0 


0 



0 

0 


/ 
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Thus in Eq. (1-11) , 



4 . Constraints 

Constraints, like objective functions, are treated the 
same regardless of whether the equations of motion are expressed 
in first or second order form. Constraints may be placed on 
state variables, derivatives of state variables, and control 
forces (accelerations, velocities, displacements and control 
forces for systems described by second order equations) . The 
general form is 

YL £ Yli + Y2u + Y3? k £ YU (1-14) 

Each row in Eq. (1-14) represents one constraint. A 
specification matrix ( MSP ) is needed to specify the nature of 
each constraint. MSP is a matrix of two columns with rows 
equal in numbers to the rows in Eq. ( 1 — 14). The elements of 
MSP are given integer values defined as follows: 
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MSP (1,1) 


MSP (1,2) 


if the i^* 1 constraint (row) is continuous 
(i.e., is to be imposed at many times) 

‘til 

if the i constraint is imposed only at 
some specific time. 


means the i^ constraints 
are of the type respectively 



Example 4 

For the same system of Example 1 suppose we want the 
constraints 


b (t) u (t) _< a (t) 

d ! 1 x 1 d 2 
y 1 d 3 

*(t 2 ) = V e 


for all t 

for t 1 only 

for all t 
for t 2 only 


In matrix form 


b (t) 

< 

— 

0 

0 

0 

0 


• 

X 

+ 

V 

fuj + 

*0 

0 


g 

< 

a(t) 

d l 


0 

0 

1 

0 


• 

y 

■ 

0 


0 

0 

. 


h (t) 

or 

d 2 



0 

0 

0 

1 


X 


0 


0 

0 



d 3 

- 


1 

0 

0 

0 


y 


0 


0 

0 



V 0 












— 

— 



- • - 
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The specification matrix is given by 

jo 0 I 


MSP = 


1 0 
0 1 


1 


3 


5 . Output 

Any linear combination of i, u, f can be tabulated 
versus time. The user inputs matrices Ql, Q2, Q3 in the 
equation ' 


Ql s + Q2 u + Q3 fy. 

The dimensions of matrices Ql, Q2_, Q3 are m x not, 
nu x not, nf x not, respectively, wKere m li” the number of 
state variables and not is the number of outputs required. 

A specific set of forcing functions, i.e., ISET(I), must be 
designated for each quantity to be tabulated. The user can 
also get trajectories of any one of the quantities tabulated. 

A trade-off curve between the objective function and 
any one constraint can be obtained by varying the bound of 
that constraint. The bounds that can be varied must be of 
the constant type (i.e., NCN = 0). For a constraint with 
both upper and lower bounds the user can either vary both 
bounds at some increment or keep one bound constant and only 
vary the other one. The details of the implementation of the 
constraint variation are given later in this Appendix. 


Example 5 

Suppose for the system in Example 2 we want to have the 
water level h. (t), h_ (t) and the controlled input e(t) to be 
tabulated as functions of time. Then 
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Q2 = 


0 

0 

1 


0 


Q3 


0 


°J., ; 

If we set ITR(l) = 1 , ITR (2) = 1, ITR(3) = 1 then the time tra- 
jectories of h 1 (t) / h 2 (t) and e(t) will be tabulated and plotted. 


Let the goal in Example 2 be to minimize: the maximum value 
of objective functions |li^(t)| and l^ft)! with the constraint 

|h^(t) | £ D. Then for each prescribed value of D we can obtain 

a minmax value of the objective functions. By varying D we can 
compute a tradeoff curve. 


6 . Input Format 

Words not underlined are variables used in the code. 
The underlined words marked with * are required in the input. 
Those underlined words marked with ** are required only if the 
data they designate are to be entered as part of the problem. 
For example, for a problem with no initial conditions (i.e. all 
initial conditions are zero), the card with INITIAL CONDITIONS 
is not needed. 

The order of the word-cards is immaterial as long as 
the appropriate data follows the proper word.. For example, the 
K MATRIX card and its data can be placed before or after the 
M MATRIX card and its data. 


7. Limitations as to Size of Problems 

The size of problems that can be solved by the 
program listed in this report is restricted by the DIMENSION 
statements in the program. The following are the current maxi- 
mum sizes of important parameters. 

NU=5, NF=5 , NSETS=1, 11=20, N0B=3, N0C=5, NOT=10, and the 
combination of the above parameters should be such that 
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Rjj < 120, C H £ 40 (See Section D of Appendix III)*. 


Card No. Columns 

Variable Format 

Comments 

1 1 

blank 


2-6 

NPB (15) 

problem number 

7-10 

blank 



11-70 

TITLE 

(12A5) 

any written title or 
description of the problem 

- / r 

71-80 

blank 

f 


2 1-5 

blank 



6 

NOE 

(ID 

st 

= 1 input system of 1 
order equations 

= 2 input system of 2 n< ^ 
order equations 

7-10 

blank 



11-15 

blank 



16-18 

blank 



19-20 

NDF 

(12) 

number of equations : of 


motion 


21-23 

blank 



24-25 

NU 

_ t , 

(12) 

number of controllers of 
the system 

26-30 

NF 

(15) 

number of elements in each set of 
forcing functions of the system 

31-35 ' 

NSETS 

(15) 

number of sets of forcing 
• functions (= 1 even if 


there are no forcing 
functions) 

* A variable dimensioned version of PREPROC is available for use 
on the NASA Langley Computer Operating System (CDC 6600) . 
Dimensional constraints are removed in an automated version of 
PERFORM at the University of Virginia Computer Science Center 
(CDC 6400) 
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Card No. 

Columns 

Variable 

Format 


36-40 

II 

(15) ' ‘ 


41-50 

TM 

(F10.6) 


51-55 

NOB 

(15) 


56-60 

NOC 

(15) 


61-64 

blank 



65 

ISP 

(ID 


66-70 

NOT 

(15) 


71-80 

blank 


3 

1-60 

ISET (t) 

I SET (2) 

• 

(6110) 


61-80 

• 

ISET (NOT) 
blank 




FORCING FUNCTION 

4 

1-16 

FORCING FUNCTION** 


17-80 

blank 


4-1 

1 

KEY 2 

(ID 


2-6 

I 

(15) 


Comments 


= 1 + number of time 
intervals for discre- 
tizing control forces 


value of time interval 

total number of objective 
functions 

total number of constraints 


= 0 minimize. the maximum 
of objective functions 

= 1 maximize the minimum 
of objective functions 

number of outputs 


where ISET (l ) designates the set 
forcing functions required by l** 1 
output ; read in up to 6 values pel 
with additional cards as needed 


= 0 continue to read 
data, e.g., cards 4-2 
are still to be read in 

= 1 last data card, e,g.,, 
if NSETS = 1, then last 
set has been reached 

the it* 1 set of forcing 
functions 


7-80 blank 
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Variable format Comments . 

KEY1 (II) ... = 0 continue to read 

data, e.g. 4-1-2 

= 1 last data card 

J (15) the j th element of i*"* 1 

set 

blank 
blank 

KEY (II) = 0 continue to read 

data 

= 1 last data card 

! 2-6 K (15) discretized time 

(1 to (II-l) ) 

7-11 . blank 

12-23 F(I,K,J) (F12.6) value of j th element of 

i th set forcing function 
at time K 

24-80 blank ' repeat until all non-zero 

elements of F(I,K,J) are 
input 

The whole deck would look 
like 

4 

4-1 4-2 

4-1-1 4-2-1 

4-1-1-1 4-2-1-1 

4-1-1- 2 4-2-1-2 

4-1- 1-3 4-2-1-3 

• • 

• • 

• • 

4-1-2 . 4-2-2 

4-1-2-1 4-2-2-1 

4-1-2- 2 4-2-2-2 
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Card No . Columns Variable Format Comments 

4-1-2-3 • 4-2-2-3 

• • , 

etc . 


INITIAL CONDITIONS 


5 

1-17 

18-80 

INITIAL 

blank 

CONDITION** 

only non-zero initial 
conditions are to be entered 

5- 1 

1 

KEY 

(ID 

= 0 continue to read data cards 
= 1 last data card 


2-6 

7-10 

I 

blank 

(15) 

the i^h element of the state 
vector § 


11-22 

23-80 

S (1,1) 
blank 

(F 12.6) 

the value of s(l,I) 
continue with card 5-2, etc. 


until all non-zero elements of 
S(1,I) are input 


5-2 has the same format as 5-1 

The initial conditions are the 
values of the response (state) 
variables at the initial time 

FIRST ORDER EQUATIONS OF MOTION (s = Al +Bu +D? k ) 

6 1-8 A MATRIX ** 

-9-80 blank 

6-1 1 KEY (II) = q continue to read data of A 

= 1 last data card of A 

•2-6 I (15) 
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Card No. Columns 


Variable Format' 


Comments 


7-11 

J 

d5) 


12-23 

A(I,J) 

(F12.6) 

value of A (I , J) 

24-80 

blank 


repeat card 6-2, etc. until all 
non-zero elements of A are input 

1-8 

B MATRIX** 


9-80 

blank 



1 

KEY 

(il) 

= 0 continue to read data of B 




= 1 last data card of B 

2-6 

I 

(15) 


7-11 

J 

(15) 


12-23 

B (I , J) 

(F12.6) 

value of B (I , J) 

24-80 

blank 


repeat card 7-2, etc., until all 
non-zero elements of B are input 

1-8 

D MATRIX** 


9-80 

blank 



1 

KEY 

(ID 

= 0 continue to read data of D 




= 1 last data card of D 

2-6 

I 

(15) 


7-11 

J 

(15) 


12-23 

D<I,J) 

(F12.6) 

value of D (I , J) 

24-80 

blank 


repeat 8-2, etc. until all non- 
zero elements of D are input. 


SECOND ORDER EQUATIONS OF MOTION (Mg + Cq + K$ + UO = Ff k ) 


9 


9-1 


1-8 M MATRIX ** 

9-80 blank 

1 KEY (II) = 0 continue to read data of M 

= 1 last data of M 
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Card No. Columns Variable Format 


Comments 


2-6 I (15) 

7-11 J (15) 

12-23 M (I , J) (F12.6) value of M(I,J) 

24-80 blank repeat until all non-zero 

elements of M are input 

10 1-8 C MATRIX ** } 

9-80 blank 


10 

1 

KEY 

(II) 

= 0 continue to read data of C 





= 1 last data card of C 


2-6 

I 

: (15) 



7-11 

J 

(15) 



12-23 

C ( I , J ) 

(F12.6) 

value of C (I f J) 


0 

CO 

1 

(N 

blank 


repeat until all non-zero 
elements of C are input 

11 

1-8 

K MATRIX* 

* 



9-80 

blank 




11-1 

12 

12-1 


1 KEY (XI) = 0 continue to read data of K 

« i 

= 1 last data card of K 


2-6 

I 

(15) 


7-11 

J 

(15) 


12-23 

. K (I , J) 

(F12.6) 

value of K(I,J) 

24-80 

blank 


repeat until all non-zero 

1-8 

F MATRIX* 

* 

elements of K are input 

9-80 

1 

blank 

KEY 

(il) 

= 0 continue to read data of F 

2-6 

I 

(15) 

= 1 last data card of F 
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Card No. 

Columns 

Variable 

Format 

Comments 



7-11 

J 

(15) 




12-23 

F(I,J) 

(F12.6) 

value of F (I , J) 



24-80 

blank 


repeat until all non-zero elements 
of F are input 

13 

1-8 

U MATRIX* 

* 




9-80 J 

blank 

- 



13-1 

1 

KEY 

(il) 

= 0 continue to read data of 

U 





= 1 last data card of U 



2-6 

I 

(15) 




7-11 

J 

(15) 

• \ 



12-23 

U(I,J) 

(Fl2.6)value of U(I,J) 



24-80 

blank 


repeat until all non-zero elements 
of U are input 



OBJECTIVE FUNCTION (PXls + PX2u + PX3fy) 


14 

1-10 

PXl MATRIX** 




11-80 

blank 




14-1 

1 

KEY 

(ID 

= 0 continue to read data of 

PXl 


= 1 last data card of PXl 
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Card No. 

Colunns 

Variable 

Format 

Comments 

• ■ 

2-6 - 

I ' • 

C(I5) 

j 


.7-11 

J 

(15) 



12-23 

PX1(I,J) 

(F12.6) 

value of PXl (I , J) . - 


24-80 

blank 


V. 

repeat until all non-zero 
elements of PXl are input 

15 

1-10 

PX2 MATRIX** 



11-80 

blank “ 

'• . 


15-1 

1 

KEY 

(il) 

= 0 continue to read data 
PX2 

= 1 last data card of PX2 


2-6 

'I 

(15) 

• 


7-11 

J 

(15) 



12-23 

PX2 (I,J) 

(F12.6) 

value of PX2 (I , J) 


24-80 

blank 


repeat until all non-zero 
elements of PX2 are input 

16 

1-10 

PX3 MATRIX** 



11-80 

blank 



16-1 

1 

KEY 

(il) 

= 0 continue to read data 
PX3 

= 1 last data of PX3 


2-6 

I 

(15) 



7-11 

J 

(15) 



12-23 

PX3 (I, J) 

(F12.6) 

value of PX3(I,J) 


24-80 

blank 


repeat until all non-zero 
elements of PX3 are input 



CONSTRAINTS 

(?L<Y_ls + 

Y2u + Y3f k <YU) 

17 

1-9 

Y1 MATRIX** 




10-80 

blank 
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Card No. Columns Variable Format 


. Comments 



1 

KEY 

(ID 

=■ 0 continue to read data 





= 1 last data card of Yl^ 


2-6 

I 

(15) 

■ . • ' ( 


7-11 

J 

(15) 



12-23 

Yl (I , J) 

(F12.6) 

value of Y1(I,J) 


24-80 

blank 


repeat until all non-zero 
elements of Yl are input 

18 

1-10 *.• 

Y2 MATRIX** 



11-80 

blank 



18-1 

1 

KEY 

(ID 

= 0 continue to read data 





= 1 last data card of Y2 


2-6 

I 

(15) 



7-11 

• J ' ' 

(15) 



12-23 

Y2 (I, J) 

(F12.6) 

value of Y2 (I , J) 


24-80 

blank 


repeat until all non-zero 
elements of Y2^ are input 

19 

1-10 

Y3 MATRIX** 



11-80 T 

blank 



19-1 

1 

KEY 

(ID 

= 0 continue to read data 





= 1 last data of Y3^ 


2-6 

I 

(15) 



7-11 

J 

(15) 



12-23 

. Y3 (I , J) 

(F12.6) 

value of Y3 (I , J) 


24-80 

' blank 

A 


repeat until all non-zero 
elements of Y3_ are input 

20 

1-21 

BOUNDS 

OF CONSTRAINTS** 


22-80 

blank 




of Yl 


of Y2 


of Y3 
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Card No . Columns 
20-1 1 

2-6 

7-9 

10 


Variable Format 
KEY1 (II) 

I (15) 


Comments - 

= 0 continue to read data card 
( 20 - 2 ) 

= 1 last data card 
the i th constraint 


blank 

MSP (1,1) (ii) = 0 ith constraint is for all 

times 

= 1 i rn constraint is for 
specific time 


11-19 blank 

20 MSP (I, 2) (II) . = 0 £ £ type constraint 

= 1 £ type constraint 

= 2 > type constraint 

= 3 _ type constraint 


21-80 blank 



MSP (1,2) 

= 

0 GOTO 

22-1AU 

■ j 

then 

go 

to 

22-1AL 

If MSP (1,1) = 0 and / 

MSP (1,2) 

= 

1 GOTO 

22-1BU 

i 

MSP (1,2) 

= 

2 GOTO 

22-1CL 


V^MSPd^) 

= 

3 GOTO 

22-1DU 


^MSP (1,2) 

■ = 

0 GOTO 

22-1EU 

' 

then 

go 

to 

22-1EL 

If MSP (1,1) = 1 and 

MSP (1,2) 

= 

1 GOTO 

22-1FU 


MSP (1,2) 

= 

2 GOTO 

22-1GL 


VJ4SP(I,2) 

= 

3 GOTO 

22-1HU 

(Use the following if MSP (1,1) = 0, MSP (1,2) 

= 

0) 



22-1AU 1 blank 
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Card No. 

Columns 

Variable 

Format 

Comments 

- ' * t 

2 

NCN 

\ 

(ID 

= 0 constant upper bound 
’ • (GO TO 22-1AU-A) 





= 1 piecewise constant upper 
bound (GO TO 22-1AU-1) 


3-80 

blank 



22-1AU-A 

1-12 

SVA1 

(F12.6) 

starting value 


13-17 

NIC1 

(15) 

number of increments (use 0 if 
there are no increments) 


18-29 

VIC1 

(F12.6) 

value of increment 

• 

30-80 

blank 


(GO TO 22-1AL) 

22-1AU-1 

1 

KEY 

(il) 

= 0 continue to read data 





= 1 last data card 


2-6 

K 

(15) 

discretized time 





NOTE: K = 2 to II if 

¥2(1,0) = 0 for all J 
K = 1 to II - 1 if 
Y1(I,J) = 0 for all J 

K = 2 to II - 1 if 
Yl (I , J) ? 0 
and 

Y2 ( I , J) ji 0 
for some J 


7-11 

blank 




12-23 

YU ( K , 1 ) 

(F12.6) 



24-80 

blank 


repeat 11-1AU-2, etc. 

22-1AL 

1 

blank 




2 

NCN 

(ID 

= 0 constant lower bound 
(GO TO 22-lAL-A) 





= 1 piecewise constant lower 
bound (GO TO 22-1AL-1) 


3-80 

blank 



22-1AL-A 

1-12 

SVA2 

(F12.6) 

starting value 


13-17 

NIC2 

(15) 

number of increments (use 0 if 
there are no increments) 




121 





Card No . Columns Variable Format Comments 

18-29 VIC2 (F12.6) value of increment 

30-80 blank 

22-1AL-1 1 KEY (II) = 0 continue to read data 

= 1 last data card 

2-6 K (15) discretized time 

NOTE: K = 2 to II if 

Y2 (I , J) = 0 for all J 

K = 1 to II-l if 
Y1(I,J) = 0 for all J 

K-2 to II-l if 
Yl(I,J) 1 0 
and 

Y2 (I , J) t 0 
for some J 



7-11 

blank 




12-23 

YL (K, I) 

(F12.6) 



24-80 

blank 


repeat 22-1AL-2, etc. 

(Use the 

following 

if MPS (I 

,1) = o 

and MPS (1,2) = 1) 

22-1BU 

1 

blank 




2 

NCN 

(ID 

= 0 constant upper bound 
(GO TO 22-1BU-A) 





= 1 piecewise constant upper 
bound (GO TO 22-lBtr-l) 


3-80 

blank 



22-1BU-A 

1-12 

SVA1 

(F12.6) 

starting value 


13-17 

NIC1 

(15) 

number of increments (use 0 
there are no increments) 


18-29 

VIC1 

(F12.6) 

value of increment 


30-80 

blank 
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Card No. Columns 

Variable 

Format 

Comments 

22-1BU-1 1 

KEY 

(ID 

= 0 continue to read data 

/ 




= 1 last data card 

2-6 

K 

(15) 

discretized time 




NOTE: K = 2 to II if 

Y2 (I , J) = 0 for all J 




K = 1 to II - 1 if 
Yl(I,J) = 0 for all J 




K = 2 to II-l if 
Y1(I,J) f 0 
and 

Y2 (I , J) ? 0 
for some J 

7-11 

blank 



12-23 

YU (K, I) 

(F12.6) 


24-80 

blank 


repeat 22-1BU-2, etc. 

(Use the following 

if MSP (I 

,1) = o. 

MSP (1/2) = 2) 

22-1CL 1 

blank 



2 

NCN 

(ID 

= 0 constant lower bound 
(GO TO 22-1CL-A) 

3-80 

blank 


= 1 piecewise constant lower 
bo.und (GO TO 22-1CL-1) 

22-1CL-A 1-12 

SVA2 

(F12.6) 

starting value 

13-17 

NIC2 

(15) 

number of increments (use 0 . 
there are no increments) 

18-29 

VIC 2 

(F12.6) 

value of increment 

30-80 

blank 



22-1CL-1 1 

KEY 

(il) 

= 0 continue to read data 




= 1 last data card 

2-6 

K 

(15) 

discretized time 




NOTE: K = 2 to II if 

Y2 (I,J) = 0 for all J 
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Card No. Columns 


Variable Format 


Comments 


K = 1 to II-l if 
Y1(I,J) = 0 for all J 

K = 2 to II-l if 
Y1 (I , J) ? 0 
and 

Y2 (I , J) ? 0 
for some J 



7-11 

blank 




12-23 

YL (K, I) 

(F12.6) 



24-80 

blank 


repeat 22-1CL-2, etc. 

(Use the 

! following 

if MSP (I 

,D = o, 

MSP (1,2) = 3) 

22-IDU 

1 

blank 




2 

NCN 

(1,1) 

= 0 constant right hand side 
value (GO TO 22-lDU-A) 


3-80 

blank 


= 1 piecewise constant right 
hand side value (GO TO 
22-1DU-1) 

22-1DU-A 

1-12 

SVA1 

(F12.6) 

starting value 


13-17 

NIC1 

(15) 

number of movements 
(use 0 if there a're 
no increments) 


18-29 

VIC1 

(F12.6) 

value of increment 


30-80 

blank 



22-1DU-1 

1 

KEY 

(11) 

= 0 continue to read data 





= 1 last data card 


2-6 

K 

(15) 

discretized time 


NOTE: K = 2 to II if 

Y2 (I,J) = 0 for all J 

K = 1 to II-l if 
Y1(I,J) = 0 for all J 

K = 2 to II-l if 
Yl(I,J) ? 0 
and 

Y2(I,J) ? 0 
for some J 
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Card No 


Columns 


Variable Format 


Comments 



7-11 

blank 




12.-23 

YU ( K , I ) 

(F12.6) 



24-80 ' 

blank 


repeat 22-1DU-2, etc. 

(Use the 

following 

if MSP (I, 

1) = 1, 

MSP (I, 2) = 0) 

22-1EU 

1 

blank 




2 

NCN 

(ID 

= 0 constant upper bound 
(GO TO 22-lEU-A) 


3-80 

blank 



22-1EU-A 

1-12 

SVAl 

(F12.6) 

starting value 


13-17 

NIC1 

(15) 

number of increments (use 
there are no increments) 


18-29 

VIC1 

(F12.6) 

value of increment 

e * « 

« h, 

22-1EL 

30r 34 
35-80 
1 

K 

' blank 
blank 

(15) 

specify the time 


2 

NCN 

(ID 

= 0 constant lower bound 
(GO TO 22-1EL-A) 


3-80 

blank 


‘ * . m ' 

22-1EL-A 

1-12 

SVA2 

(F12.6) 

starting value 


13-17 

NIC2 

(15) 

number of increments (use 
there are no increments) 


18-29 

VIC2 

(F12.6) 

value of increment 


30-34 

K 

(15) 

specify the time 


35-80 

blank 



(Use the 

following 

if MSP (I, 

1) = 1,MSP(I,2) = 1) 

22-1FU 

1 

blank 




2 

NCN 

(ID 

= 0 constant upper bound 
(GO TO 22-lFU-A) 


3-80 

blank 



* 

22-1FU-A 

1-12 

SVAl 

(F12.6) 

starting value 

- 

13-17 - 

NIC1 

(15) 

number of increments (use 
there are no increments) 


18-29 

VIC1 

(F12.6) 

value of increment 
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Card No 


Columns 


Variable Format 


Comments 



30-34 

35-80 

. K { I ) 
blank 

(15) specify the time 

(Use the 

following 

if MSP (1,1) 

= 1, MSP (1/2) = 2) 

22-1GL 

1 

blank 



2 

NCN 

(11) = 0 constant lower bound 

(GO TO 22-1GL-A) 


3-80 

blank 

* < 

22-1GL-A 

1-12 

SVA2 

/ 

(F12.6) starting value 


13-17 

NIC2 

(15) number of increments (use 0 if 
there are no increments) 


18-29 

VIC2 

(F12.6) value of increment 


30-34 

K 

(15) specify the time 


35-80 

blank 

" '■ 

(Use the 

following 

if MSP (I ,L) 

= 1,MSP (1,2) = 3) 

22-1HU 

1 

blank 



2 

3-80 

NCN 

blank 

(11) = 0 constant right-hand side 

(GO TO 22-1HU-A) 

22-1HU-A 

1-12 

SVA1 

(F12.6) starting value ' 


13-17 

NIC1 

(15) number of increments (use 0 if 
there are no increments) 


18-29 

VIC1 

(F12.6) value of increment 


30-34 

K 

(15) 


35-80 

blank 

OUTPUT (QlS + Q2u + Q3f k ) 

23 

1-9 

10-80 

Ql MATRIX** 
blank 
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Card No. 

Columns 

Variable 

Format 

Comments 

23-1 

1 

KEY 

(ID 

= 0 continue to read data 





= 1 last data card 


2-6 

I 

(15) 



7-11 

J 

(15) 



12-23 

Q1(I/J) 


value of Ql (I / J) 


24-80 

blank 


repeat until all non-zero 
elements of Ql. are input 

24 

1-9 

Q2 MATRIX** 



10-80 

blank 



24-1 

1 

KEY 

(il) 

= 0 continue to read data 





= 1 last data card 


2-6 

I 

(15) 



7-11 

J 

(15) 



12-23 

Q2(I,J) 

(F12.6) 

value of Q2 (I , J) 


24-80 

blank 


repeat until all non-zero 
elements of Q2 are input 

25 

1-9 

Q3 MATRIX** 



10-80 

blank 



25-1.. 

1 

KEY 

(il) 

= 0 continue to read data 





= 1 last data card 


2-6 

I 

(15) 1 



7-11 

J 

(15) 



12-23 

Q3 (I , J) 

(F12.6) 

value of Q3 (I , J) 


24-80 

blank 


repeat until all non-zero 
elements of Q3^ are input 

26 

1-10 

TRAJECTORY** 



11-80 

blank 
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Card No. 

Columns 

Variable 

Format 

Comments ' ; 

26-1 

1 

KEY 

(ID 

= 0 continue to read data 
e.g. 26-2 





= 1 last data card 


2-6 

I. 

(15) 

the i^* 1 output 


7-10 

blank 




ii 

ITR 

(il) 

= 0 no plot wanted 




i 

= 1 want to plot iD* output 


12-80 

blank 

* .. / 

repeat until all needed 
trajectories are specified 

27 

1-17 

TRADE OFF 

DIAGRAM** 


18-80 

blank 



27-1 

1-4 

blank 




5 

ITD 

(ID 

= 0 no plot wanted 

' 




= 1 want tradeoff diagram 


6-9 

blank 




10 

KSWT 

(ID 

= 0 the ordinate of tradeoff 


curve starts from 0 
= 1 the ordinate of tradeoff 


curve starts from the minimum 
value of objective function 

11-80 blank 


VERIFICATION OF INPUT AND CALCULATED,. DATA 

. ( » , ' • . 


28 

1-6 

VERIFY** 


7-80 

blank 

i — i 

00 

CN 

1-4 

blank’ 


5 

INP 


*1 * 


(ID 


= 0 does not print input data 
= 1 print out input data 


128 


Card No. Columns Variable Format Comments 



6-9 

blank 



10 

INT 

(11) = 0 does not print data 
calculated by PREPROC 


11-80 

blank 

= 1 print data calculated by 
PREPROC 

■ J 


FINISH INPUT DATA 

29 

1-17 

FINISH 

INPUT DATA* 


18- 8 0 1 

blank 

MODIFICATION 

30 

1-14 

MODIFY 

PROBLEM** 


15-80 

blank 

Put the modified data between 
cards 30 and 31. The user can 
modify all but the data in the 
first three cards. The format 
of the modified data are the 
same as described before. 

31 

1-18 

FINISH 

MODIFY DATA** 


19-80 blank 

■.■I • 


RESTARTING A NEW PROBLEM 


32 1-7 RESTART ** 

8-80 blank ' Put data from card 1 through 

card 31 after this card. This 
permits the user to run two. or 
more different problems at a 
time. 


END OF THE PROBLEM 


33 1-4 STOP * Use this card when there is 

no more data to be read 


5-80 blank 
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APPENDIX II 


LINEAR PROGRAMMING FORMULATION OF THE 
LIMITING PERFORMANCE PROBLEM 


In this appendix the limiting performance problem is 
formulated as a linear programming (LP) problem. We consider 
dynamical systems in which those portions for which the per- 
formance is to be measured are represented by time dependent 
control forces . Although the control forces can replace non- • 
linear system segments, the remaining portions of the system 
must be linear. In addition, the kinematics of the whple 
system must be linear. Although the engineer frequently 
describes a dynamical system with second order differential 
equations, it is more convenient for us to work with first 
order equations. Hence, Section A of this appendix considers 
the conversion from second to first order equations. 

The unknown control forces are discretized in time. 
Integration of the equations of motion for known disturbances 
yields the response (state) variables as linear combinations of 
the control forces. The objective functions and constraints 
comprising the limiting performance problem are formed from 
these response variables and hence are also linear functions 
of the discretized control forces. The limiting performance 
problem then can be formulated as a linear 1 programming problem. 
Section B of this appendix treats the integration of the first 
order equations of motion containing the unknown control 
forces. Finally, we formulate the limiting performance prob- 
lem as a linear programming problem in Section C. 
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- A. Conversion of Equations of Motion of a Dynamic System to a 
System of First Order Equations 

Consider the general equations of motion of a multi-degree-of 
freedom system in the form 


+ C<5 + Kq + Uu = Ff (II-l) 

where 

q = displacement vector (n x 1) 
u = control force vector {nu x 1) 

I = forcing function vector (nf x 1) 

M = mass matrix (n x n) 

C = damping matrix (n x n) 

K = stiffness matrix (n x n) 

U = matrix of constants associated with control forces 
(n x nu) 

F = matrix of constants associated with external disturbances 
(n x nf) 

Letters with bars over them refer to vectors, e.g., 


q = 


^1 

< 12 . 

3n J 


is the displacement vector. 


while letters with bars under them refer to matrices, e.g.. 


C = 


C 11 

c 12 * 

c 

rH 

o 

• 

• 

• 

C 21 

• 

• 

CM 

CM 

o 

••• c 2n 

• 

fnl 


* * • c nn 

equations 

can be 

of 

first 

order « 


is the damping coefficient matrix. 
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2nxnf 

ft 

Then the system of second order equations can be written as 
Ws+Xs+Yu=Zf 

(II-2) 

• 

or s=As+Bu+Df 



where 
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0, 0^ and 0 2 are null matrices of order nxn, nxnu and nxnf 
respectively. 

][ is an nxn identity matrix. -Here we have been made use of 


W 


-1 _ 


■M" 1 C 


M 


-1 


M" 1 


i 

i 

. 1 . 


M"1 


i 

i 0 

i — 


-J 2nx2n 


Note that A, B, D, W and W - ^- are partitioned matrices. 


B. Integration of a System of First Order Equations 


We wish to integrate the system of first order equations 


s = As + Bu + Df { 1 1 — 3 ) 

with initial conditions 

_l q(o) 

s < t - 0 = 5 (o) = [ ] 

q(o) 

The vectors s, u, f, and matrices A, B, D are defined in the 
previous appendix . 

To solve Eqs (II-3) let 

i (t) = e~ t z ( t ) . ( 1 1 — 4 ) 

Substitute Eq. (II-4) into Eqs. ( II— 3 ) to find 
. dz 

Ae- z + e- 11 — = Ae^ 1 z + Bu + Df (II- 5) 

- dt - “ 

If we substract Ae^-2 from both sides of Eq (II-5) and then 

, -1 

premultiply both sides by e'^'we obtain 

_ e "At [Bu + Df] 
dt 
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and z(0) = s(o), since s (o) = e^. lo, z(o) = lz(o) = z(o). 
Then, 

2 ( t ) = s(o) + e _ — ^ [Bu (y) + Df(y)]dy (II-6) 

We see from Eqs. (II-4) and (II-6) that the solution of the 
system (II-3) can be written as 

S(t) = e^tKo) + S t e“£y [Bu (y) + Df(y)ldy] (II-7) 

O *" 


where y is a dummy variable for integration. 

To complete the integration in Eq. (II-7) we must know u(t) 
and f(t). The loading functions f(t) are prescribed. The 
forces u(t) are unknowns to be found as solutions of an 
optimization problem. We discretize u(t) so that between two 
discretized times u is represented by constants. This_will 
permit Eq. (II-7) to be integrated. _ For convenience, f(t) is 
discretized in the same fashion as u(t). 


Piecewise Constant Discretization 


Assume u is discretized as a piecewise constant function 
of time as shown in Fig. II-l. Here, between t = iT and 
t = (i + 1 ) T , u(t) = u(iT), where T is the time increment At. 


Then, at any time kT 

S(kT) = e£. kT [i(o) + / kT e~- y [Bu (y) + Df(y)]dy] 

k-1 (i+1) T 

= e& kT s(o) + e^ T I f e-Ay[Bu(y) + Df(y)]dy 

i=o iT 


Now, 

e"^y[Bu(y) + Df(y)]dy 


= /.(i + D T e - ^y [Bu (iT) + Df(iT)Jdy 
IT — — 

= e -&(i+l)T A -l[e& T - I] [Bu ( iT) + Df(iT)] 

Thus , 

s(kT) = e^ kT s (o) 

+ e A ( k ~ i~ 1 ) T A - 1 ( e AT - I)[Bu(iT) + Df(iT)] (II-8) 

i=o 
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is the required discrete formal integration of Eq. (II-3) > for 
a piecewise constant representation of u(t). Note, here 

A - 1 ( e^ T - 1) = A -1 (1. + AT + - — + . + ••• - I) 

“ 21 3! 

= IT + (II-9) 

21 3! 4! 

Even if A is singular, the solution, i.e. , Eq. (II— 8) , 'is 
still valid since the series on the right-hand side of Eq. 

(I I- 9) always exists. ' 

C. Formulation of the Linear Programming Problem 

In this Section the limiting performance problem is formu- 
lated as a linear programming (LP) problem, with discretized 
control forces as unknowns. This is accomplished using the 
objective functions and constraints in the form described in 
Appendix I together with the solution of the equations of 
motion given in Section B of this Appendix. 

The solution of the first order equations of motion are 
given in Section B for the piecewise constant discretization 
of the control forces as: 


s(kT) = e^ kT I (o)tl 1 eA( k -i-D T A _1 (e^ T -I) [B u(iT) + D f (iT) ] 

1-0 ( 11 - 10 ) 

In order to be compatible with indices in F0RTRAN , let i(l) be the 
initial condition S (0) . We rewrite equation (11-10) as: 

k-1 

3 (k) = (e^Tjk-lgd) + | (e£ T ) k - 1 - i A- 1 (e£ 1, - I)[Bu(i) + D f (i) ] 

(11-11) 

where s(l) = 3(0) 

s(k) = s( (k-l)T) .. 

Now, define the following matrices: 
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AE I = A -1 (e^ T -I) 
R(n) = P n_1 AEI B 
T (n) = pH" 1 AEI D 


( 11 - 12 ) 


Note that R ( 1 ) = AEI B, T(l) = AEI D 

Here R(i) and T(i) (i = l,2,...n) are matrices of order m x nu, 
m x nf, respectively, where m = number of first order equations 
to be integrated and n is the number of discretized time intervals. 
In terms of those new matrices Eq. (11-11) becomes 

s(k) = P k “ 1 s(l) +*Z*[R(k-i)u(i) + T (k-i) f (i) ] (II-13) 


Now, define 

RG (k) = [R(k-l) IR (k-2 ) :•••: R(l)] 

TG (k) = [T (k-1) : T (k-2) • • • • j T (1) ) 
"" r ™ • 


f (k) = 


?U) 

hv 


f (k-l) 



u (k) 


u(l) 
• • • 
u (2) 


Q (k-l) 


Using Eqs. (11-14), Eq. (11-13) .becomes 


I (k) = P k_1 I(l) + RG (k) U (k) + TG (k) f (k) 


(11-14) 


(11-15) 
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Let 

C (k) = cT(k) + C2(k) 
where cT(k) = P^sd) 

C2 (k) = TG(k)f(k) 

We can simplify Eq. (11-15) to 

I (k) = C(k) + RG(k)u(k) (11-18) 


For a system subject to several sets of forcing functions, 
Eq. (11-18) becomes 

s(s,k) = C(s,k) + RG (k) u (k) (11-19) 

where 

C(s,k) = CT (k) + cl(s,k) 

Cl(k) is given in Eq (11-17) 

C2(s,k) = TG(k)f (s,k) 


f(s,k) = 


Eq. (11-19) expresses the state variable vector s as a 
linear function of the control forces, initial conditions, and 
forcing functions. Here the initial conditions and forcing 
functions are prescribed, the only unknowns are the control 
forces (u) . 

Since the objective functions and constraints are all 
formed as linear combinations of s, u, f, they are all linear 
functions of the u's. Then the problem of finding the u's that 
minimize the maximum or maximize the minimum value of the ob- 
jective functions subject to constraints becomes a problem in 
linear programming. We will now define the limiting performance 



' (11-16) 
(11-17) 
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problem in terms of the LP problem and explain how to set up 
the data in the proper form. 


The LP problem is to find z such that c z is minimized (or 
maximized) subject to HzfG (11-20) 

Where c = (1, 0, 0, ,0) is a row vector 


z = 


with u (II) = 


u (I I], 

§ 111 . 

u(2) 


(H-21) 


ji<II-l)J 

H and G are a matrix and a vector constructed from the 
objective functions and constraints, respectively. 


The objective function is formed as 
? = PXls + PX2u + PX3f (s) 


( 11 - 22 ) 


where the index s designates the s th set of forcing functions 
are being used. There are two cases of interest: minimizing 

the maximum value of $ or maximizing the minimum value. The 
minimization of the maximum value of 9 is equivalent to finding 
the minimum of <p such that 


| V± | £ <J> for all i = 1 to NOB 

or - <J> + <_ 0 (11-23) 

<J> + Vi > 0 (11-24) 

where is the i fc ^ row of V 

Using Eqs. (11-19), (11-22), Eqs. (11-23) , (11-24) become 



Eqs. (11-25), (11-26) are valid for k = II, II + 1, ... L 

s = 1,2, NSETS 

where PXlj., PX2^ , PX3_ i are the i th rows of the matrices PXl , 

PX2 , PX3 -respectively . The values of II and L depend on whether 
PXl-j or PX2^ are zero. If PXl, = (0 0 ... 0) then II = 1, 

L = II-l. If PX2i= (0 0 ... 0), then 11=2, L = II. If 
PXl^ (0 0 . ..“57 and PX2^ * (0 0 ... 0), then II = 2, L = .11-1. 
The contribution to H and G for a typical objective function 
(e.g., i*-* 1 row of y) can be found from Eqs. (11-25) and (11-26). 
It is put_in matrix form in Table II-l. The total contribution 
to H and G from the objective function is a matrix containing 
NOB submatrices in the form of Table (II-l) , where NOB is the 
number of objective functions. 

To maximize the minimum value of ¥ is equivalent to 
finding the maximum of 4> such thatl^i^ <{>. Note that to find 
maximum is equivalent to finding minimum -<p. Thus the LP 
problem for this case is still in the form of Eq. (11-20) except 
now the row vector C is (-1,0,0, .... 0). The formation of 
arrays H and G for the objective functions remain the same as 
above except the inequality signs are reversed. 


The constraints are formed as 
Yl S + Y2u + Y3f (s) 

The type of i th constraint depends on the value of the i^ row 
of the specification matrix, i.e. MSP(i,l) and MSP(i,2). The 
specification matrix has been defined in Appendix I. For 
MSP(i,l) = 0, the constraint is applied at all times. 


1. MSP ( i , 2 ) =0 

The constraint is 

YL i <Yl i i + Y^u + Y3J? (s^YU^ (11-27) 


140 


. 5 ( 1 ) 
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Using Eq. (11-19), Eq. (11-27) can be written as ‘ 

Yl i RG(k)5(k) + Y2iu(k)<YU i (k) - Yl^s-fk)- - Y3.jf (s>k) 

(11-28) 

YliRGfkJufk) + Y^u (k) ^YL^k) - YO^Cfs, k) - Y3jf (s,k) 
for K = XI to L 


The contribution to H and G of Eqs. (11-28), (11-29) are summarized 
in Table II-2A, and Table II-2B, respectively. The value of II and 
L depends on whether YljOr Y2. is zero. If Yl- = (0 0 ... 0), 
then II = 1, L = II-l. If Y2 j = (0 0 ...0), then 11=2, L = II. 

If Ylj ^ 7 * (0 0 ... 0) and Y2-; ? (0 0 ... 0), then II = 2, L = II-l. 


2. MSP (i , 2) =1 

The i'" 4 * * * * * * 11 constraint is 

\ Yl-j 8 + Y2.J.U + Y3 i f(s)<YU i (11-30) 

The contribution to H and G is the same as Eq. (11-28). 

3. MSP (i , 2) = 2 

Here the constraint is expressed as 

Yli§ + Y^U + Y3 i £(s)>YL i (11-31) 

The contribution, to H and G is the same as Eq. (11-29) . 


4. MSP (i , 2) = 3 

In this case the constraint is given by 

YljS + Y2 a u + Ya^fts) = (11-32) 


The contribution to H and G is the same as Eq. (11-28) , except 

now the inequality becomes an equality. For MSP(i,l) = 1/ the 

constraint is applied at a specified time only. The formal 

contribution to H and G for all- cases is the same as for 

MSP(i,l) = 0, except k is now restricted to a specific time, 
i.6. , KS ( l ) . 
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TABLE I I- 2 A 





D. Dimension of the Linear Programming Problem 

The dimension of the linear programming problem depends 
on the number of controllers of the system, the number, and 
nature of objective functions and constraints, the number of 
time intervals, and the number of sets of forcing functions. 

It is independent of the number of degrees of freedom of the 
system. The dimensions of the matrix H and vector G can be 
calculated as following: 


R h = [ (II-l) x (2N X + 2C 1 + C 3 ) + (II-2) x (2N 2 + 2C 2 + C 4 > + 

2C 5 + Cg] x NSETS 

C H = 1 + Nu x (II-l) 

R G = R H 

where 

R h = Number of rows of matrix H 
C H = Number of columns of matrix H 
r g = Number of rows of vecter G 

N^ = Number of objective functions for which pxi-j = 0 or PX2^= 0 , 

N 2 = Number of objective functions for which PXl ^ ^ 0 and PX2 ^ ^ 0 

C-^ = Number of constraints for which 

MSP (1,1) = 0, MSP (1,2) = 0 and = 0 or Y2 ± = 0 
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= Number of constraints for which 

MSP (1,1) = 0, MSP (1,2) = 0 and Yl^ ^ 0 and Y2± ? 0 

C 3 = Number of constraints for which 

MSP (1,1) = 0, MSP (1-2) t 0 and Y2 i = 0 or .= 0 

for which 

t 0 and Yl ± t 0 and Y2 ± ? 0 

= Number of constraints with MSP (1,1) = 1,MSP(I,2) » 0 

Cg = Number of constraints with MSP(I,1) ,= 1, -MSP (1,2) ^ 0 

II = Number of time intervals used (see Appendix I) 

NSETS = Number of sets of forcing functions 

As an example, consider the dimensions of the linear programming 
problem for the SDF problem shown on page 17 . Here 11-10, N 2 = 1, 

Ci = 1, NSETS = 1, Nu = 1 

R H = [ (10-1) x (2 x 1 + 2 x 1) ] x 1 = 36 

C H = 1 + 1 x (10-1) = 10. 


C 4 = Number of constraints 
MSP(i,l) = 0 ,MSP (1,2) 
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APPENDIX III 


PROGRAMMING DESCRIPTION OF THE COMPONENTS OF PERFORM 


This appendix contains descriptions of the three major 
components Of PERFORM: PREPROC, MPS/360, PSTPROC. 


A. PREPROC 


This section includes a description of program PREPROC. 
Appendix IV A contains a listing of program PREPROC. 

Table III-l contains a listing of the arrays used in the 
program with explanatory notes. 

The variables beginning with an IVE card followed by a 
number are used to communicate to the printing routine those 
arrays which contain nonzero elements or special conditions. 
Table III-2 relates IVE numbers to program arrays. 

1. Main Section - 

The main section of the program is described in the 
PREPROC section of Section III of this report. 

2. Subroutines (Alphabetical Order) 

a. BODATA 

This subroutine punches the BOUNDS action of the 
MPS/360 linear programming problem. The default bounds for 
MPS/360 variables are zero lower and unlimited upper; because 
solutions to limiting performance problems can be negative, 
the default overriding BOUNDS section is necessary. 

b. BRIDGE 

_ This subroutine calculates the. arrays R, T, cT, 

C2, and C as described in Appendix II. It is called from the 
main section. 


4 
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Array 

Identifiers 

A 

B 

CAY 

CC 

CM 

CYM 

D 

EM 

EMI 

F 

FF 
FM 

G 
H 


Dimension •' 
m* x m 


m* x nu 

n x n 
n x n 
n x n 

n x n 
m* x nf 

n x n 
n x n 


Notes Common Blo< 


First order system A is 
input. Second order A 
computed by subroutine 
COVERT 

First order B is input. 
Second order computed 
by subroutine COVERT 

Matrix K in model 

Matrix C in model 

Becomes M~^C in sub- 
routine COVERT. 

Becomes M“^K in sub- 
routine COVERT. 

First order D is input. 
Second order computed by 
COVERT . 


'■V ABD 


ABD 

EKC 

EKC 

INV 

INV 

ABD 


Matrix M in model EKC 

Becomes M“1 in . subroutine INV 

MIV. 


nsets x II x nf Subscript 1: forcing 

function set number, 2: 
time interval , 3 : row of 

equations of motion with 
which associated. This is 
not model matrix F. 


n x nf 


Matrix F in model. 


n x nf 


Becomes M"^F in subroutine 
COVERT 


number of LP 
rows (See 
Dimension of 
. LP problem) 


Right-hand side of LP 
tableau 


number of LP Linear programming 
rows; number of tableau. 

LP columns (See 
Dimension of LP 
problem) 


FOFUN 


EKC 

EKC 


RHV 


PVAL 


. TABLE III-l., PREPROC ARRAYS 

♦ 
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Array 

Identifiers 

Dimension 

Notes Common Block 

I SET 

not 

■ r ISET (1) identifies forcing 
function of the set to be 
used in tabulation row I 

TAB 

ITITL 


Alphanumeric problem title 

TIT 

ITR 

not 

ITR (I) is 0 if plot of I th 
output is not desired. 
Otherwise plot is made 


IXX 

noc 

IXX (I) indicates the behavior 
of the bounds of 1^ constraint 
= 0 no varying bounds 




= 1 YU varies 




= 2 YL varies 




= 3 YL and YU vary 


KS 

noc 

'KS(I) specifies the specific 
time for I th constraints 


MSP 

noc x 2 

Specific type of constraint 
See Appendix I for description. 

ISP 

NRE 

number of 
LP rows 

Identifies relation of each row 
with its right-hand side 

RHV 



"G" means >_• 


PX1 

nob x m 

"L" means <_ 
"E" means = 
* * 

i 

OBJ 

PX2 

nob x nu 

* * 

OBJ 

PX3 

nob x nf 

* * 

OBJ 

Q1 

not x m 

* * * 

TAB 

Q2 

not x nu 

* * * 

TAB 

Q3 

not x nf 

■ * * * 

TAB 

S 

2n 

Used for initial condition 
in program 

INIT 

UM 

n x nu 

Becomes M~^U in subroutine 
COVERT 

INV 

UU 

n x n 

Matrix U in model 

EKC 


TABLE III-l. PREPROC Arrays (continued) . 
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Array 

Identifiers Dimension 


Notes 


Common Bloi 


Y 

Y1 

Y2 

Y3 

YU 


YL 


C 


P 


member of Array of varying bounds , 

varying bounds PSTPROC tradeoff diagram 

noc x m * * * * 


, r . UPLO 
CONST 


noc x nu 


* * * * 


CONST 


noc x nf 


* * * * 


CONST 


II x noc 


Upper bound vector. UPLO 

Subscript 1: time interval 

2 : number of constraints 


II x noc Lower bound vector, sub- UPLO 

scripts have same meaning 
as above. 

nsets x II x m Arrays pertinent to the forma- 
tion of right-hand side of the 

. .. .. .. < . LP problem. - (See Appendix II). 

Subscript 1: forcing function CCS 

set 

2: time 

3: row number in state 

(response) vector 

m x m Becomes e^L in BRIDGE 


C2 


Arrays used in form array C. 


Cl II x m 


Arrays used to form array C. 


R II x m x m Arrays pertinent to expressing BRG 

responses as function of the 
control forces. 


T 


AEI 


II x m x nf 


m x m 


Arrays pertinent to expressing BRG 
responses as functions of 
forcing functions. 

Becomes -A -1 (e£ fc - I)in BRIDGE MYAP 


TABLE III-l. PREPROC Arrays (continued) 
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Symbols used in array table . 

nob = Number of objective" functions (rows in PXl , PX2 , PX3 
expression ). 

noc = Number of constraint expressions. 

nsets =Number of sets of forcing functions. 

not = Number of tabulation expressions. 

nf = Number of forcing functions. 

nu = Number of isolators (control forces). 

* m = NDF For first order system, m = 2 (NDF) for second 
order system. 

** Arrays related to input of objective functions, 
see Appendix I . 

*** Arrays related to input of required output, 
see Appendix I . ' 

'*** Arrays related to input of constraints, see Appendix I. 

NDF Number of equations in the system of equations. 

TABLE III-l. PREPROC Arrays (concluded) 


i 
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IVE 


Array Read or Command 


1 

2 

3 

4 

5 

6 

7 

8 
9 


. 10 
11 
12 

13 

14 

15 

16 


17 


18 

19 

20 
21 


A 

B 

C 

EM 

CL 

CAY 

FF 

UU 

S (initial condition) 

F (forcing function) 

PX1 

PX2 

PX3 

Y1 , 

Y2 

Y3 

MSP 

(Not used) 

Q1 

Q2 

Q3 


TABLE III-2. IVE Variable Array Relation 
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c. 


CODATA 


This section punches the COLUMNS section of the 
MPS/360 file. The coefficients of the tableau array H are in 
this section, CODATA is called by LPDATA which is called by the 
main section. 

d. COVERT 

This subroutine reduces the equations of a 
second order system to those of a first order system. It is 
called from the main section. 

e. ELEMTP 

ELEMTP is called from the main_section to compute 
the objective function portions of the H and G arrays. 

f . ELEMTQ ' : 


The constraint portions of the H and G arrays 
are computed by subroutines called by ELEMTQ which is, called 
from the main section. The routines called depend on whether 
the programming tableau rows are equalities or inequalities 
and whether YL, YU, or both are part of problem constraints. 

g. ELEMTR 


Subroutine ELEMTR computes the appropriate 
portions of the H and' G array for the problem constraint I 
which is a passed parameter from ELEMTQ. Table II-2A of 
Appendix II shows the computation of ELEMTR. As can be noted, 
this routine computes only the YU array. 

h. ELEMTS 

_ This subroutine completes portions of the H 

and G arrays dependent on problem constraint I which is a passed 
parameter of ELEMTQ. Table II-2B of Appendix II shows the 
computation of ELEMTS. Note that this routine works only with 
YL. 


i. . ELEMTT 

Routine ELEMTT is called when the constraint of 
passed parameter I is an equality. Note that equalities always 
involve the YU vector. 
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LPDATA 


j • 


The subroutine LPDATA punches section header 
cards for the MPS/360 linear programming problem and calls other 
subroutines to punch the actual section records. Subroutines 
RODATA, RHDATA, BODATA, and CODATA are called. LPDATA is 
called from the main section. 

k. M12X3 


This subroutine multiplies a two-dimensional 
array times a portion of a three-dimensional one, both of which 
are parameters . 

l. M2 2X2 

This routine multiplies two passed arrays. 

m. MEXAP 

This subroutine evaluates the matrices P and 
AEI used in Section C of Appendix III. These matrix functions 
are computed using approximate formulas described in Ref. 9. 
This routine can be replaced by any routine that uses other 
algorithms to compute matrix exponential functions.” 

n. MIV 


This routine computes the inverse of array EM 
and stores it in array EMI. 


o. POST 

Subroutine POST punches the PSTPROC report 
spec ifica tio ns . It includes arrays Ql, Q2, Q3, Y, R, F, C, 
ISET and ITC. No subroutines are called by POST. 

p. PRNT 

After the linear programming problem_has been 
com piled subroutine PRN T is called to print the H, G, and 
NRE arrays. Array NRE shows the sign (<>=>>) of the linear 
programming row. 

q. RHDATA 


Subroutine RHDATA is called by LPDATA to punch 
the RHS section of the MPS/360 linear programming problem. 
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r . RODATA 


Subroutine RODATA punches the ROWS section of 
the MPS/360 linear programming problem. The ROWS section 
associates names with the tableau rows and specifies the re- 
lation of each row to its right-hand side, i.e., , = , or >. 

S. SEARCH 

This subroutine sets a parameter to 1 if there 
are nonzero elements in the given row of the passed two- 
dimensional array. Otherwise a zero is returned. 

t. VERFY 


If a VERIFY indicator record is encountered and 
variable INP is not zero, subroutine VERFY is called to print 
all of the PERFORM problem specification data. 

u. VERFY 2 

If variable INT is nonzero, subroutine VERFY2 
is called to print out the intermediate resulting arrays. 

These include A, B, D, R, AEI . 


B. MPS/360 

A listing and discussion of the MPS/360 control program 
is in this section. The material about MPS and its relation 
to PERFORM in part B of the Section III of this report is useful 
background information to this appendix. 

Consider the listing of the MPS/360 control program. 

The first two statements, PROGRAM and TITLE, identify the 
program name and the title to be printed on each page of the 
solution report. 

The first instruction is the macro INITIALZ. Macro 
instructions direct the compiler to include at tie point of 
appearance a predefined set of source instructions identified 
by the macro name; macros eliminate the need for tedious 
programmer coding of repeatedly used sets of instructions. 

The INITIALZ macros initializes values of certain tolerances, 
standard frequencies, etc, and specifies actions to be taken 
should errors occur during execution. 
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MPS/360 Control Program Listing 


//UVAAE001 JOB 605R0.WANG ,MSGLEVElxi 


/♦FORMAT Pll >nnNAMExSYSPUNCH»DESTxLOCAL 
/♦MAIN T IMF * 1 O.LINFSxl 0 
// EXEC LP 
//lp.sysin on * 

PROGRAM TIPMII 

T I T L r t» UVA EXAMPLE5 ll 

INITIAL/ 

MOVE r XOAT A# IEXMP1 I ] 

MOVE rXPBNAMF» ILTMOYNI ] 

CONVERT t » CHECK I >* SUMMARY I) 

8COOUT tlONEH 

SETUP f IMINI , IROUNO * * * BDS 1 I ) 

PICTURE 

MOVE rxORJ, IOBJFI I 
MOVE r XRHS» IRHS1 » 1 
PR IMAI 
SOLUTION 

MOVE T XOATA. I ANSWER* 1 
PUNCH t *L I ST I # I VALUE * I 
EXIT 
PENO 

/* 

//GO.SYSPRINT DO SYSOUT * A 
//GO.SYSIN 00 * 

(Place MPS/360 Linear Programming Problem Data Deck here) 
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The next two move instructions set data cells XDATA and 
XPBNAME. Identifiers beginning with X are predefined special 
variables for passing information to the procedures. The con- 
tents of XDATA, which is set to EXMP1, identifies the nam-> of 
the data to be used by the following procedures. If the input 
stream has several named decks one can be picked for use by 
setting XDATA to the name on the deck NAME card. XPBNAME is 
the name to be given to the problem file — an internal MPS 
file, omitted from the system flowchart and designated PROBFILE. 

Procedure CONVERT writes on PROBFILE in converted for- 
mat the data deck named by XDATA. The contents of PBNAME 
identifies the converted data. Parameters 'CHECK' and 'SUMMARY' 
specify special checks to be performed and a summary of the data 
to be printed . , * 

The next procedure, BCDOUT, prints the problem on the 
PROB file identified by XPBNAME. Parameter 'ONE' specifies 
one value per printed line. . 

Certain computations that must be done before the solu- 
tion is found are done by SETUP; XPBNAME identifies the problem 
to be set up. Parameter 'MIN' specifies object function mini- 
mization; 'BOUND' and 1 BDS ' identifies" BDSI as the vector name 
of the bounds specifications. 

PICTURE causes a picture of the matrix to be printed. 

The next two move statements identify vector’ OBJF to be 
used as the objective function coefficients and RHSl the right- 
hand side vector. 

PRIMAL is the optimization procedure. The next procedure 
SOLUTION produces a printout of the solution obtained by pro- 
cedure PRIMAL. , 

The linear programming problem solution card file is 
produced by procedure PUNCH. Parameter 'LIST' produces a 
printed listing of the output and 'VALUE' is necessary for 
values to be punched. 

Procedure EXIT is the final statement and returns con- 
trol to the system from the program. It is equivalent to the 
FORTRAN CALL EXIT or STOP statements. 

PEND, like FORTRAN END, specifies to the compiler the 
end of the program. 
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C. PSTPROC 


A discussion of the program PSTPROC is given in this 
section. Appendix IV B contains a listing of program PSTPROC. 
There are three subroutines in PSTPROC in addition to the main 
section. Each of these is described below. 

1 . Main Section 

. The main section is described in the Section III of 
this report. Arrays have the same identifiers in both this 
program and PREPROC. 

2. Subroutines 


a. SEARCH 

Subroutine SEARCH sets a parameter to 1 if a 
nonzero element is found in a given row of the passed array. 

If all row elements are zero, a 0 is returned. 

b. TRAJ 

Subroutine TRAJ plots the tabulation function 
against time; a tabulated listing is printed in addition to the 
line printer plot. 

The logic of the program is straightforward; 
the position of the star is determined and computed GO TO 
statements transfer control to a WRITE -FORMAT pair which prints 
the star. 


c . TRADE 

The logic of this subroutine is much like that 
of TRAJ. It plots the tradeoff between the objective function 
values and the values of varying bounds. 
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APPENDIX IV 


LISTINGS OF PROGRAMS 
A . Program PREPROC 
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000163 GUP) = 0. 
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